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Abstract 

A systematic evaluation method of matrix elements of quantum one-body op- 
erators between nucleon states in the three-quark Bethe-Salpeter(3qBS)/Faddeev 
approach in the NJL model is reviewed. We do not confine ourselves to a partic- 
ular truncation scheme of the 3qBS kernel. One of our main aims is to derive a 
general condition to be imposed on a given BS kernel in order that the PCAC rela- 
tion is satisfied correctly. We apply this condition to some particular 3qBS kernels. 
We numerically calculate g-j^NN in the 3qBS /Faddeev approach to estimate the 
violation of the Goldberger-Treiman/PCAC relation due to the UV-regularization 
scheme. We finally consider the non-vanishing current quark mass effects on the 
Goldberger-Treiman relation. 
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1 Introduction 



The Nambu-Jona-Lasinio(NJL) model is a quantum field theoretical effective quark model 
based on QCD. While the confinement, which is one of the most important properties 
of the low energy QCD, is not incorporated, it is the simplest quantum field theoretical 
quark model respecting the chiral symmetry and providing us with explicit examples how 
the spontaneous breaking of the chiral symmetry plays an important role in hadronic 
phenomena m, in particular, the interpretation of pion, K, and rj as Nambu-Goldstone 
modes. Based on these advantages, not only mesons but also baryons have been exten- 
sively studied in the NJL modelp], |^. The studies of baryons are mainly classified into 
the following two categories: (CI) the relativistic mean field approach[Q, and (C2) the 
relativistic three-quark Bethe-Salpeter(3qBS )/Faddeev approach]^, ||, |^, |^. Since the 
NJL model is a second quantized field theory, it is, in principle, possible to consider the 
effects of qq excitations in the nucleon. This is one of the most interesting targets in 
the studies of the nucleon structure beyond the non-relativistic constituent quark models. 
In this respect, on the one hand, the meanfield approach incorporates qq effects through 
the mesonic hedgehog fields leading to the solitonic picture of the nucleon, where, how- 
ever, the manifest Lorentz covariance is unfortunately missing. On the other hand, the 
3qBS /Faddeev approach, which respects the manifest Lorentz covariance leading to the 
quark-diquark picture of baryons[0], approximates the 3qBS kernel according to the lad- 
der truncation scheme, which prevents us from studying non-trivial qq effects beyond the 
excitations of the "three-quark RPA vacuum" ||TT|(the backward diagrams included in the 
ladder approximation). 

It is thus necessary to extend the 3qBS kernel in order to study the non-trivial qq 
effects in a Lorentz covariant manner, and therefore attempts have been (and are) made 
to go beyond the ladder truncation scheme |jl2|. To extend the NJL 3qBS kernels, chiral 
symmetry imposes important constraints. In particular, it would be very useful to have 
a criterion which tells us which kinds of kernels lead to the PCAC relation^ correctly. 
One of the main aims of this paper is to derive such a criterion and to provide explicit 
examples of its application to some of the existing 3qBS kernels. We will first review 
how matrix elements of quantum one-body operators are evaluated systematically in the 
3qBS framework by introducing classical external fields as a technical tool. Note that 
the evaluation of the bound state matrix elements in the 3qBS framework is not obvious 
beyond the diagramatic argument. We will provide a direct formula to evaluate the 
matrix element in terms of the Faddeev framework. We then consider which properties 
of the 3qBS kernel are required to satisfy the PCAC relation correctly. We consider the 
relevant Feynman diagrams to calculate the matrix elements for several 3qBS kernels, and 
apply our criterion to these 3qBS kernels. Although we are going to restrict our attention 
only to the PCAC case, these considerations themselves can be extended to other cases 
straightforwardly such as the electromagnetic current, the isospin current, the baryon 
number current, etc — actually, the PCAC case is the most complicated one0. 

^We mean by "PCAC relation" in the sense of Eq.(^). 

^The proof of the electromagnetic Ward identity in the Faddeev approach to the NJL model in the 
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Even if the truncation scheme is consistent with the PCAC relation, the UV regular- 
ization schemes, which have been adopted so far in practical numerical calculations in the 
NJL model, usually, spoil it. It is thus worth while to evaluate Qt^nn and qa in the chiral 
limit in order to estimate explicitly the violation of the Goldberger-Treiman relation due 
to the UV regularization scheme. We will see that with the Euclidean sharp cut-off the 
violation is up to 4 %. The extracted value of g-^NN is 13.2. We will also estimate the 
PCAC violation due to the cutoff scheme off the chiral limit, which is again up to 4%. We 
will evaluate the off-shell gT,]^^{= q-^nn) with the single-pole dominance approximation 
for the form factor hA^q^)- We obtain a very reasonable value q-^nn = 13.5. However qa is 
6 % larger than the experimental one. We will also consider the effect of the non-vanishing 
current quark mass on the Goldberger-Treiman relation, and find that it is in the opposite 
direction than expected from the experimental values. We will analyze this by including 
also the effect of the qq interaction in the color singlet iso-vector axial-vector channel, 
providing an analytic expression for the deviation from the Goldberger-Treiman relation 
by assuming that the vacuum is approximated according to the mean-field approximation. 
We will find that the effect of this additional qq channel works into the unwanted direc- 
tion. We will see that, to resolve this problem, it is necessary either to improve the gap 
equation for the vacuum beyond the mean field approximation or to evaluate the on-shell 

Q-kNN- 

2 Bound State Matrix Elements and the Criterion 
(Sufficient Condition) for the PCAC Relation 

We consider effective quark Lagrangians with global SU{2)f x 5'f/(3)c symmetries: 

C = i^{i^-m^)ij + Ci, (1) 

where mo is the current quark mass, and tp and if) are the quark bispinor fields. Ci is 
a local, chirally symmetric four-fermionic interaction Lagrangian of NJL type, i.e., Ci = 

gr (^ipTtjj^ , where F is a matrix with Dirac, isospin, and color indices. Examples are 



r 

\2 tTT.^, ;^„u\2 



the original NJL type Cj = g \^{^p^p) — (ip'y^Tijj) j and the color current interaction type 

^1 = 9 ("^7^^^) 5 etc. (For detail, see Appendix 0.1.) Since interaction Lagrangians 
of this type can all be dealt with in the same manner, we will not confine ourselves to a 
particular one. 

To evaluate a matrix element in the 3qBS framework is not really straight forward. 
This is because what is directly obtained in the 3qBS framework is not the nucleon eigen- 
ket |A^) but the 3qBS amplitude^, which itself is a matrix element of the type {OlTipipiplN) . 

ladder truncation scheme can be found in ref.j^], which was apphed to the quark-diquark model in 
psf . Although these authors did not use the external field method, the results are consistent with our 
formalism. A systematic approach using the external field method can be found in ref . . 

■^The 3qBS amplitude is often referred to as the "wave function" due to historical reasons. However, 
we prefer to call it as "3qBS amplitude" to avoid unnecessary confusions pq] . 
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We might thus suspect that a huge amount of information, which is contained in the ket 
vector |A^), could be missing in the 3qBS amphtude. In the relativistic Faddeev frame- 
work, the situation is shghtly more comphcated — there are no immediate representations 
of the Faddeev amphtude in terms of the canonical operator formalism^. Therefore, we 
will first review how to evaluate a matrix element in the 3qBS framework based on the 
canonical operator formalism. 

To this end, it is convenient to introduce space-time dependent external fields: 

= vl{x) — , a^{x) = a'^^ix) — , m{x) = s{x) + ip''{x)-t^ — , (2) 

where v^{x) is an isovector vector field, a^(x) is an isovector axial- vector field, s{x) is 
a scalar isoscalar field, and p"(a;) is a pseudo-scalar isovector field. (The iso-spin index 
a runs over 1, 2, 3.) The following Lagrangian density characterizes how these external 
fields couple to the associated quantized one-body operators: 

£W = + A (3) 

where V^^ = tjj'j^—ip and A"''^ = tpl'^J^—i'- The superscripts [e] or [0] will be used to 

indicate quantities in the presence or absence, respectively, of the external fields v,a,m — 
rriQ. For simplicity, we include the current quark mass mo in the definition of the external 
field m. We assume that these external fields, v,a,m, are localized in space-time, i.e., 
they vanish except for some finite space-time region. After performing the necessary 
manipulations, we will eventually set a{x) = v{x) = and m = mo in the whole space- 
time. This limit will be refereed to as the "vanishing external fields limit" . 

The 3qBS equation in the presence of the external fields is derived from the Schwinger- 
Dyson(SD) equation[0 by means of the same technique as for vanishing external fields: 



= + G!fVWG["], (4) 

where V^t^' is the sum of 2PI and 3PI interactions, and Gq^ is an anti-symmetric combi- 
nation of products of three constituent quark propagators in the presence of the external 
fields. For simplicity, we denote 

To avoid cumbersome notations, we often adopt the operator notation to suppress the ex- 
plicit integration symbol J (i'^xi 6/^X2(^^X3, the explicit space-times coordinates {xi,X2,X3) 

and the explicit indices of the quark fields. 

Since the external fields are localized in space-time, a nucleon state "propagates" 
freely in the "past". Therefore, G^^l has an asymptotic initial nucleon pole. The residue 

"^The Faddeev amplitude provides the same amount of informations as the 3qBS amphtude does, as 
we show in Appendix 
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at this nucleoli pole satisfies the following homogeneous 3qBS equation in the presence of 
the external fields: 

where \l/[Ai"^(P)]['^] is the 3qBS amplitude for the nucleon in the presence of the external 
fields with the isospin a, helicity a and asymptotic four momentum P = {Ei\i{P'^), P) 

with E]\}{P'^) = \l m\i + P2 {iTiN '■ nucleon mass). In terms of the canonical operator 
formalism, it is expressed as: 

nK{P)f\x,,x,,x,) (7) 







X exp -i / d^x vl{x)V''^'{x) + <(x)A"'^(a;) + ^ (m(x) - mo) ^ K{P) 



0) and 



N^{P)j are the vacuum and the nucleon state vectors, respectively. They are 

eigenstates of the "unperturbed" Hamiltonian if'"' obtained by setting y = a = 0, m = 
ruQ. We adopt the covariant normalizations: 







O) = 1, lNi,{P') N^{P)) = M^!!) (2vr)35(3)(p' _ P)S^,J^,^ 



rriN 



(8) 



Note that quantized operators are represented in the "interaction picture" in Eq.(^, where 
the "unperturbed Hamiltonian"^ is if'^l and the "perturbing Hamiltonian" is 



Hr = 



d^x (v1^{x)Vaf,{x) + a)^(x)A''^(x) + V(x) (m(x) - mo) ^(x)) . 



(9) 



We next consider the matrix element (A^|y4|^|A^). Applying the functional derivative 
iS 

to both sides of Eq.(P), we obtain the following relation in the vanishing external 



Sa^j^ix 



field limit, 

'iS-^[N;^{P)P 



5a^(x) 



5a^(x) 



(10) 



The expressions for the "3qBS amplitudes" entering in Eq.(|TO]) are given in terms of the 
canonical operator formalism as follows: 



■^[N^{P)f\x,,X2,Xs) 

^N^{P)f\x^,X2,X3) 

l5^[N^{P)P{^U^2,X3) 







Tijix,)ij{x2)ij{xs) K{P) 



fill 



iV-(P) Ti,{x^)ij{x2)^P{x^) 







T^(xi)7/>(a;2)V'(x3)^'^(a:) K{P) 



Note that H^^'^ is the fuh Heisenberg Hamiltonian in the absence of the external fields. 
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By using the canonical operator analysis, i.e., by inserting into the Green's function the 
complete set relation in the baryon number 1 sector: 



d^p rriN 



K{p))(K{p) 



+ 



;i2) 



and by using the Fourier transform of the step function 9{t), we see thatQ the product 
\I/[o]vl>'''' provides the residue of the Green's function at the nucleon pole, i.e., 

d'^p rriN i 



{2'nYEN{f)po-EN{p')+ie 



(13) 



In Appendix ^, we explain how to obtain and normalize these two 3qBS amplitudes based 
on the relativistic Faddeev framework. In Appendix 0, we prove that the residues '^^^^ 
and satisfy the following 3qBS equation in the absence of the external fields(cf Eq.(||)) 

irM^[°l[iV^(P)] = ^[°l[iV^(P)], t™[iV;j(P)]ir[°l = t'°'[Ar^(P)], (14) 

where a tilde is used to indicate an "amputation", i.e., *'[iV;j(p)][°l = '^[N^{p)]^^'^G^Q^''^ . 
By rearranging Eq.(|l^), we obtain the following relation: 



1 - KM vH(^) 



On the one hand, near the nucleon pole the resolvent is given as 
1 



1 - i^M 



^[0]^[0]-l 

f d^p rriM 



(15) 



(16) 



{2TiYEN{f)p''-EN{p')+te 
The l.h.s. of Eq.(^), on the other hand, is given near the nucleon pole as follows: 



(17) 



E 

a' a' 

E 



d^p rriN 



{2TiYEN{f)p^-EN{f)+te 
d'^p rriN i 



i2n)^ E^ip") p"" - E^if) + te 
By comparing both sides of Eq.(^), we are left with: 







MM 



^{Nl{P)f\ 



^The argument is similar to the one given in p. 92 in Ipl 



A similar consideration leads to the following relation: 



^V'(a;)75y^(a;) 




A combination of these two relations provides us with: 



(19) 



(20) 



^k;(p' 




— 2mi 



i5 



Sp'' 



XI 



These expressions show that the 3qBS amplitudes provide enough informations to evaluate 
these matrix elements. Here the "relative time" dependence of the 3qBS amplitude surely 
plays a very important role. At any rate, we should keep in mind that, although bound 
state matrix elements are obtained by sandwiching one-body operators in the canonical 
operator formalism, it is these two-body operators shown on the r.h.s. of Eq.(^D|) that 
should be sandwiched in the 3qBS framework'', which leads to crucial differences ||14|. 
The corresponding formula to evaluate the bound state matrix elements in the Faddeev 
framework is derived in Appendix ^ Note that it is not hard to extend these arguments 
beyond the NJL model. 

If the full 3qBS kernel were at our disposal, the chiral symmetry of the original La- 
grangian should directly lead to the PCAC relation in the 3qBS framework. However, 
because the 3qBS kernel is in practice truncated such as to be manageable, it may hap- 
pen that the truncation scheme spoils the chiral symmetry. We are thus interested in 
the criterion to decide which 3qBS kernel gives rise to the PCAC relation correctly. For 
this purpose, it is convenient to introduce the local "axial" gauge transformation of the 
external fields as follows: 

z^-/")(x) -^(")(x)75 -m(^) = Q{x){i^-]^{x)-^{x)-f5-m)Q{x), (21) 



n{x) 



2 ' 



where and m^^^ are the gauge images of f , a and m, respectively. We assume 

that uj{x) is also localized in space-time. Infinitesimally, these gauge transformations are 
expressed as follows: 



6u!'^{x) 
6u!'^{x) 
Suj'^ix) 



,,,a'(y)5^'\x ~ y) 



(22) 



dx^^ 



5('\x - y) + eatcv'Jy)5^'\x - y) 



lp^{y)6^'\x~y) 



^If 3PI interactions are included in the 3qBS kernel, there appear three-body operators in addition to 
the two-body operators in the 3qBS framework. 
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where eabc is the totally anti-symmetric tensor. Let F[f,a,m] be a functional. We define 
the functional derivative - — through the following relation: 



5F _ 5F[t;(^),a('^),m(^)] 



(23) 



By using the chain rule, 

OU!'^[Xj 

tives with respect to the external fields: 



i6 



6uj^{x) 



is expressed as a linear combination of functional deriva- 

(24) 



+ 



^ Ssjy) _^ 6p\y) iS 



i6 



6iu^{x) Sa^^{y) 6iu^{x) Sv'^{y) 5uj''{x) 6s{y) 6uj^{x) Sp'^{y) 

-2s{x 



i5 



In particular, we have the following relation in the vanishing external field limit: 



5uj^{x) 



— 2m 







^5p^{x)^ 



(25) 



where [e^'^^] denotes the gauge transformed fields Now we state the suf- 

ficient condition for the validity of the PCAC relation as follows: 
Sufficient Condition: 

If a 3qBS kernel behaves in the following manner under a gauge transformation f2(x): 

i^['''"'l(xi,X2,X3;yi, 7/2,2/3) (26) 

then the PCAC relation is satisfied. 

Our "criterion" is obtained from Eq.(P^ by applying the functional derivative with 
respect to the gauge transformation on both sides, and using Eg. (p5D . We are left with 
the following chiral Ward identity, which is the "criterion" for the validity of the PCAC 
relation: 

Criterion: (The Chiral Ward Identity for the 3qBS Kernel) 



i6 



''Salix 



— 2m, 



i5 



6p^{x) 



(a;i,a;2,X3; 2/1,2/2, 2/3) 

[(rS5<S>l®l)5^^^(a; - 2/1) + (l(8)r''75(g)l)(5(^)(x - 2/2) + {l0l®T^l5)S^'^\x - y^) 

{T%(g)l0l)6^^\x - Xi) + {10T%01)6^^\X - X2) + il0l(g)r^To)S^^\x - X3 



(27) 
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To see that this indeed leads to the PCAC relation, we sandwich Eq.(^) between 
'^[N^',{P')f^ and ^N^{P)p and use Eq.(|l4D and Eq.(0) to get 

{n^,{P') \d^A'^{x) - 2monP{x)'y5T'^{x)\N^{P)) = 0, (28) 

which is the desired relation. 

We should comment here on the reason why we called Eq. ( p6D a "sufficient condition" . 
If we try to extend our method to SU{3) l x SU (3)r-QCD in the presence of external chiral 
gauge fields, the relation (^) is no longer valid due to the existence of the non-Abelian 
anomaly^ [|l6l. Note, however, that, even if there is such a non-Abelian anomaly, because 
the anomalous contributions are polynomials of at least second order in the external chiral 
gauge fields, the infinitesimal form (the criterion Eg . (|27|) ) is still valid. However, if more 
than two functional derivatives of the external chiral gauge fields are involved, one must 
pay full attention to the non-Abelian anomaly even for the infinitesimal form. At any 
rate, the chiral SU{2)l x SU{2)r case, which is of most interest to us here, is known to 
be anomaly free. 



3 Feynman Diagrams 

The aim of this section is to see which kinds of Feynman diagrams are involved in the 
expression Eq.(|T^) of {N\A^^{x)\N) , before applying the criterion/sufficient condition to 
particular 3qBS kernels. 

We ffist consider the constituent quark propagator S^p{x^y) in the presence of the 
external fields. The self energy is, in principle, obtained from the sum of IPI diagrams. 
However, in practice, it is approximated in the mean field (Hartree-Fock) treatment, as 
expressed by the following gap equation: 

iSf{x,y) = tS§A^,y) (29) 

- / d^ziSl^^F{x,z)USP{z,y)TT{TtSPiz,z)) . 

r 

gp^^ are the effective coupling constants in the qq channels, which are defined in Appendix 
0.1. sl)}p{x,y) is the current quark propagator, which is defined through the following 
relation: 

(z^(^) -^(x)75 -^(a;) - mix)) ^-Ua^, 2/) = ' v)- (30) 

The solution to Eq. (|29D is obtained as a self-consistent solution to the following equations: 

S^ix^y) = (^^-^75-^-m-SW)"'(a;,y) (31) 

sW(z) = ^2#Tr(4^1(z,^)r)r. 
r 

^This anomaly should not be confused with the anomaly, i.e., the Abelian anomaly. 
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T}-^\z) is the self-energy of the constituent quark. Now, several comments are in order. 
(1) Due to the presence of the external fields, there is no translational symmetry any 
more, and therefore the self-energy depends on the space-time point z. (The dependence 
on just a single space-time coordinate z is due to the mean field approximation.) (2) Non- 
vanishing external fields may lead to not only a non- vanishing scalar condensate, but also 
a non- vanishing pseudoscalar, vector, axial vector condensates. However, in the vanishing 
external field limit, we have mo + S^"^ = M (the constituent quark mass). (3) Because 
the NJL model is non-renormalizable, it is necessary to introduce a UV regularization in 
all loop integrals. Hereafter, whenever a loop integral is encountered, the integration is 
understood clS cL TG gularized one. (4) If this regularization respects the chiral symmetry, 
S^p{x,y) and T}-^\z) should behave in the following way: 



S[^'^'l(z) = Q{z)T.^''\z)n{z). 



(32) 



For later convenience, the transformation of Sp is depicted in Fig.|l|. 

All diagrammatically truncated 3qBS kernels K'^^^ can be expressed as the product 
of several constituent quark propagators S^p and several elementary local vertices F. 
Because, in our case, these elementary local vertices do not depend on the external field^, 
K^^^ depends on f , a and m only through S^p. In a symbolic notation, this may be denoted 
as i^t^l = /^[F; S"!;^], and the functional derivatives are then symbolically expressed by 
using the chain rule as follows: 



J2 J dhd^z' 



aa' 



Z, Z 



6K[T-Sf] 



(33) 



where a, a' are triples of the Dirac, iso-spin, and color indices. This relation implies 
that the general rule to obtain the functional derivative of the 3qBS kernel is to replace 
each constituent quark propagator in turn by its derivative, and then sum up the resulting 
terms. Now all we need is 6S^p / 6a!'^{x), which is obtained by applying 5/5a^(x) to Eq.(PTD 
as follows: 



5ai{z) 



5T^^\z' 



+ / d^z'Sf{x,z') 



sP{z,y) 



5SW(z') 



r 



/ 


'ssP{z',z') 


[0] \ 




5al{z) 




V 







Sf{z',y). 



F. 



(34) 



(35) 



^If derivative couplings are involved in the interaction Lagrangian, our argument is straightforwardly 
extended by replacing ordinary derivatives by chiral covariant derivatives. In this case, 5/6a''^{x) also hits 

the argument F in if [fI^^I; 5*15;']. This is essential to extend the arguments given in the next section. In 
this case, it is not simply A^f^{x) but the conserved current operator, that is inserted in the r.h.s. in the 
3rd line in Eg. (pi]) — A''^^{x) fails to be a conserved current in this case. 
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By inserting the first relation into the second, we obtain the following closed equation for 
6J:/5al. 



Y: 2#Tr ( z) [ 1 S^H^. z')V ] V 



(36) 



+ j d^z"Y.2tgf'^i:i [5[°l(z',z") 



5Y.^^\z"\ 



5M(z",z')r r. 



For simplicity, we adoptj^ the effective qq coupling constant in the iso-vector axialvector 
channel (y^ax = 0. (For the precise meaning of (^ax, see Appendix 0. We will discuss the 
general case, i.e., ^^ax 7^ 0, in Appendix |^.) We can parameterize our solution as follows: 



/ 



d^z'e"^^''-'^ 



5T.^%z') 



(37) 



and from Eq. (|36|) we see that H{q^) satisfies the following equation: 

H{q^) = -2tg^IiUq^)-2ig^Ii,^{q^)H{q^), 

where Ii^A{.<f) and 1155(5^) are the bubble integrals which are defined in Appendix 
The solution is expressed geometric series: 



(3^ 



11.2. 



H{q' 



-2ig^Ii^A{q^) 
l + 2^(?^n55(g2)' 



(39) 



5Sp / 5a!'^{z) is depicted in Fig.|^.(a). The two diagrams on the r.h.s. correspond to the 
two terms in Eq . (pi]) , respectively, where the second one is proportional to q^H{q^) and 
has the pionic pole. Note that, whenever the quark propagator is obtained as a non- 
trivial solution to a self-consistent equation which leads to the chiral symmetry breaking, 
SSp/Sa^^ always has a non-trivial "mesonic part" H{q^). 

Now we can consider particular 3qBS kernels and the associated diagrams relevant for 
the matrix element calculation. Our first example is the ladder truncated 3qBS kernel, 
which is depicted in Fig.^(a). Solid lines represent the constituent quark propagator. A 
slash indicates amputation of the constituent quark propagator. The interaction strengths 
in the various qq channels are obtained by applying the Fierz identity to the interaction 
Lagrangian. (See Appendix |B[) In this case, Eq . (|33|) leads to the diagram in Fig.||(b), 
where the insertion on the quark line indicated by the "®" has been defined in Fig.|^. 
(Diagrams with the same topologies are omitted.) 

The second example is the qq interaction involving the exchange of qq pairs in the 
t-channel ("meson exchange interaction"), which is depicted in Fig.^(a). In this case, 
Eq.(P5D leads to the diagrams depicted in Fig.^(b). We see that, in addition to the 



""^^The non-vanishing effective qq coupling constant in the iso-vector axial-vector channel gax does not 
change the chiral symmetry properties [pJ. 
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coupling of the external field to the external quark propagators, couplings to internal 
quark propagators are also involved (the second diagram), which are often referred to as 
the "meson exchange current" contributions. 

It is straightforward to extend these considerations to more complicated and realistic 
cases. A moral is that, once the 3qBS kernel is specified, a unique set of Feynman diagrams 
exists to determine the matrix element of the axial vector current. 

4 Applications of the Criterion/Sufficient Condition 

The aims of this section are to see whether the 3qBS kernels presented in the previous 
section satisfy the criterion Eq. (^) /sufficient condition Eq. ([26|) , and to provide examples 
how to use the general considerations of Section 2 practically. 

We first consider the chiral transformation properties of the elementary local ver- 
tices in the Lagrangian. The global chiral symmetry of the interaction Lagrangian 

Lj = '^gri'ip'^ipY implies the following identity: 
r 

J2g^ (e^^-^e^re^^-^e^)^ ^^.....e^e^^.^e^)^^ = J^g.T^^T,, (40) 
P *j P 

Since Cj is a contact interaction without any derivative terms, this identity remains valid 
even if Q'' acquires space-time dependence, i.e., 

J2gr {n{x)m{x))^, {n{x)m{x))^^ = Y^grT^.Tki. (41) 
P u P 

Actually, these qq interactions consist of several "closed chiral multiplet sectors" , and the 
identity of the type Eq.(pD holds in each such sector separately. For example, the qq 
interaction in the vr and a mesonic channel forms a closed chiral multiplet under the local 
axial gauge transformation, i.e., the following identity holds: 

S^Ai-j:{l,T^),,{l,T^),i = {n{xMx))^,{n{xMx))^^ (42) 

a=l 3 

- ^ {n{x)j,T^n{x)),, {n{x)j,T^n{x))^^. 

a=l 

We need to establish similar relations in the qq sector. For this purpose, it is convenient 
to use the following Fierz identity (for details, see Appendix |B[l.): 

= Eari^r = E#^(^r'v^^)(^'^rV). (43) 

r r' 

In this representation, the identity Eq.(^) is expressed in the following way: 

{nix)mxf) {nixfv'nix)) = E#^r:/',, (44) 
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These qq interactions also consist of several closed chiral multiplet sectors, and in each 
sector Eg . (|4^) is valid separately. In particular, the qq interaction in the scalar diquark 
channel {J'^ = 0~^, isoscalar, color 3) forms a closed chiral singlet, i.e., 

E {{l5C-')T2f3A)^.{{C^,)T,f3A)^^ (45) 



A=2,5,7 

= J2 (^(^) ((75C-^)r2/3^) n{xf).,{nixf {ic^,)T,/3A) nix] 

A=2,5,7 



kl 



where is the rescaled color Gell-Mann matrix (3a = \J\^a with the normalization 
ti{j3A(3B) = 35a_b- Note that Pa for A = 2,5,7 are anti-symmetric matrices corresponding 
to the color 3c diquark channels. The qq interaction in the axial vector diquark {J^ = l"*", 
isovector, color 3) together with the vector diquark {J'" = 1~, isoscalar, color 3) channel 
forms a closed chiral multiplet, i.e.. 



E 

A=2,5,7 



E {il,C-'){T^T,)(3A) {{Cr){r2ra)(3A)^^ (46) 

. a=l 

- {il,l6C-')r2PA)^^ ((^757^)^2/9, 

^ {Q{x) {{i,C-'){rar2)(3A) n{xf)..{n{xf {{CY){r2ra)(3A) ^{x))^ 

. a=l 

- {n{x) {h,7,c-')T2(3A) n{xf)^,{n{xf ((^757^)^2/3^) ^(x))^ 



E 

A=2,5,7 



We give a list of these closed chiral multiplets in Appendix 0.1. 

Now we consider whether the sufficient condition Eq.(PBD is satisfied in the case of 
the ladder truncation scheme, i.e., the 3qBS kernel of the type Fig.^ In this case, the 
2PI interaction consists of only an elementary contact interaction. In practical numerical 
calculations, it is further truncated according to the quantum numbers of the diquark 
channels, i.e., the scalar diquark channel, the axialvector diquark channel, etc. We assume 
that the truncated vertex corresponds to a sum of closed chiral multiplets. Fig.^ shows the 
steps involved in the analysis of the gauge transformation properties of the kernel. In the 
first step, we apply the local gauge transformation Q{x) in order to get the transformed 
kernel (l.h.s.of Eq.(^)). Only the constituent quark propagators transform, and their 
transformation is given by Eq.(^2|). In the second step, we use Eq. (Pf) , and the last 
step is due to the fact that amputated propagators are delta functions. As a result, the 
condition Eq.(^) is satisfied by the kernel in the ladder approximation, provided that the 
truncation of qq channels is done such as to have closed chiral multiplets. In particular, 
since the qq interaction in the scalar diquark channel forms a chiral singlet, the 3qBS 
framework in the ladder truncation scheme keeping only the qq interaction in the scalar 
diquark channel gives rise the PCAC relation correctly. However, if the qq interaction in 
the axialvector diquark channel is further included, it is in principle necessary to include 
also the qq interaction in the vector diquark channel. In practice, however, the vector 
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diquark channel can be expected to have small effects, at least on the nucleon mass, from 
the nonrelativistic analogy. 

Next we consider whether the sufficient condition is satisfied in the case of the kernel 
involving qq exchange, i.e., the 3qBS kernel of the type Fig.^ Fig.^ shows the argu- 
ment for one of the infinite terms involved in the ladder sum. We assume that each 
vertex corresponds to a sum of closed chiral multiplets. In the first step we apply the 
gauge transformation using Eq.(|52D to transform the constituent quark propagator. In 
the second step we use Eq.(^). Note that all the phase factors, which appear around the 
internal vertices, cancel themselves. The last step is due to the fact that the amputated 
propagator is just a delta function. This demonstrates the validity of Eq.(^) for the 
kernel involving qq exchange. Several comments are in order. First, since the mass of the 
pion is so small, the qq exchange in the pionic channel will contribute much more than 
the one in the sigma mesonic channel. However, in principle, it is necessary to include 
both channels in order to get the PCAC relation correctly. Second, if we include the 
qq exchange interaction in the 3qBS kernel, it is necessary to take into account also the 
"meson exchange current contributions" in the calculation of the matrix element of the 
axial current, since the l.h.s. of our criterion Eq.(p7D involves also the second diagram in 
Fig.^.(b), as discussed in Section 3. 

These examples are straightforwardly extended to more general, complicated and re- 
alistic cases (for example, to the expansion scheme of We may consider these 3qBS 
kernel as a formal sum of 2PI and 3PI diagrams. It is thus not hard to convince ourselves 
that the chiral symmetry in the original interaction Lagrangian leads to the PCAC rela- 
tion correctly in our formalism, provided the truncation of the kernel is done consistent 
with the condition (^). 

5 The Goldberger-Treiman Relation and Qt^nn 

In this section, we give the explicit numerical results for Qt^nn in the chiral limit together 
with the off-shell Qt^nn off the chiral limit. We restrict our attention to the ladder trun- 
cation scheme keeping only the qq interaction in the scalar diquark channel. Although 
this truncation scheme of 3qBS kernel should lead to the correct Goldberger-Treiman 
relation, the UV regularization (the Euclidean cutoff in our case) unfortunately violates 
the chiral symmetry. One of the aims of this section is therefore to estimate the degree 
of violation of the PCAC/Goldberger-Treiman relation due to the UV- regularization. In 
this section, unless explicitly indicated for the more general case, the external fields are 
understood to he absent, i.e., f = a = 0, m = mg. For consistency reason, we prefer to 
use the axial current operator A^^{x) as an interpolating field for the pion in this paper. 

For those readers who prefer to use ip{x)'~)^—ip{x) as an interpolating field for the pion, 
the following relation would be convenient|^: 

K^) = f.d.Trlix) + ■■■, (47) 

^^We will not use this relation explicitly in this paper. 
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where 7r^g(x) is the asymptotic field operator for the pion. 



5.1 The Extraction of Qt^nn in the Chiral Limit 

In this subsection, tuq = is understood. The axial form factors gA{<f) and hA{q^) are 
defined through the following relation: 

^{")a(p)^ + q^hA{q^)l,) u^^^L) = (iV(")'^(P) \a'^{x = 0) | iV(^)^(L)) , 

(48) 

where u^"^°'(P) and u^^^'^{L) are the final and initial nucleon bispinors with the isospins 
a, c, the helicities a, [3 and the momenta P, L (P^ = = m^), respectively, q = P — L 
is the momentum transfer, qa and QnNN are extracted as follows [p!5|: 



9A = 9Aiq' = 0), M9') = ^%^ + ---- (49) 

q^ 

A multiplication of q^ on both sides of Eq. (^S|) leads to 

= ^(")'^(p)^(^^(g2)^^^ + g2^^(^2)^^)^(/3)c^^) 

= u(-)-(P)I^ {2m^g^{q')j, + q'h^iq')^,) u^^^^L) 

= 2m^gA{q^)+q^hA{q^). (51) 

where we used the PCAC relation (^) for itlq = 0. As is well-known ||15||, due to the pion 
pole of h^iq^) at q^ = 0, this relation leads to the Goldberger-Treiman(GT) relation: 

rriNgA = firg-KNN- (52) 

We expand the both sides of Eq.(^8D with respect to a small q^ to extract Qt^nn- The 
l.h.s. is expanded as follows: 



^(")a(p)L. ^g^(q^)Yj, + q^hAiq'H) u^^HA^'P) (53) 



ui-)-{pf- [gA{q')l^l, + q'hA{q')l,) S{A~')u^P^^{P) 

r-b \ / ^b 



where (A^^)^*" = {6-^)^," is a boost matrix, i.e., = {A-^)^''P^, S{A''^) = e^^'^^ is 
the boost matrix in the Dirac bispinor space with S{—X) = {—X)^^—a^'^. Note that 
q^ = p^- {A-^)^''P^ - X^^P^ + 0(A2). We thus have A = 0(g). Hence we have 

(«^"^"(^)y75^(-A)«(^)^(P)) q^ = Oiq'), (54) 
15 



which cancels the pion pole of hA{q^) at = 0. 

The r.h.s. of Eq.(^) is evaluated by the 3qBS /Faddeev expression Eq.(|n 
Eq.(p3D of the nucleon matrix element: 



and 



jV(°)"(p) A^^'ix = 0) iV(^)'=(L) 



5ai{x = 0) 



i5K 



5a»(x = 0) 



miiL)] (55) 



where tp[NI^{L)], ip[N^{P)] and Ky are the Faddeev amplitudes and the Faddeev kernel, 
respectively. These quantities are defined in Appendix ^, where also their relation to the 
3qBS quantities is explained. The Faddeev equation is depicted in Fig.|^. Since we adopted 
the ladder truncation scheme in this section, the qq interaction in Fig.|^ is understood to 
be point-like and separable. Hence the double line in Fig.0 can be expressed as a geometric 
series of qq bubble integrals. In addition, since the qq interaction is truncated to the scalar 
diquark channel, the double line in Fig.^ now stands for the following t-matrix tsd('?^ 
the scalar diquark channel: 



m 



tsd(g') = 2 X 



2ig, 



sd 



1 - 2igsdnsd{q'^ 



with 



nsd(g')5 



AB 



tr {{C'y5r2PA)^SF{q + k){^,C-\2l3B)iSF{-kf 



(56) 



(57) 



(2vr)4- 

Since the second line in Eq.(^) involves the functional derivative of the Faddeev kernel, 
the contributions are classified as the following three types0(cf. Appendix |A.3|) : (a) the 
quark current contribution[Fig.||.(a)], (b) the exchange current contribution[Fig.|[(b)], (c) 
the diquark current contribution[Fig.||.(c)]. Note that, due to the iso-scalar nature of the 
scalar diquark, the diquark current contribution to the matrix element of the iso-vector 
axial current vanishes identically. 

The evaluation of the diagram Fig.^.(a) posses a problem, since there appears a delta 
function 5^^\p — I — q/2) in addition to the delta function associated with the total mo- 
mentum conservation. Therefore, depending on the momentum transfer q, the evaluation 
involves the values of the Faddeev amplitude at points which are outside the mesh used for 
the solution of the Faddeev equation^. To avoid this problem, we use the homogeneous 
Faddeev equation Eq.(^)[Fig.|^ to iterate the Faddeev amplitude in the final state. The 
diagram Fig.^(a) is thus exactly equivalent to Fig.|^.(a'), which is free from the additional 
delta function. We therefore have to evaluate the diagrams of Fig.H.(a') and Fig.|[(b). 
Since the operator insertion on a constituent quark line has been given in Eqs. (B3) - 



^^In the diagrams of Fig. 8, we adopted the relative momenta of the spectator quark and the diquark, 
which are defined in Eq.(121), with the value rj = 1/2. 



^■^The numerical procedure to obtain the mass and the associated Faddeev amplitudes is explained in 
detail in |, 0, |. 
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(p9|) and Fig.0.(a), we are led to the four diagrams shown in Fig.||.(d). Here, we used the 
identity Eg .(11271) to express the matrix element only in terms of the Faddeev amplitudes 
in the rest frame, and used also the identity 

= S{A)^S{A-'). (58) 

Note that, since we truncated the qq interaction to the scalar diquark channel, two of the 
boost matrices ^(A"^) out of the three (5'(A^^)®'^) in Eq. ( |127| ) cancel. The remaining 
S{A~^) is found in Fig.H.(d) at the point where the operator is inserted, which is a 
consequence of the manipulation Eq. (|58|) . By using the following relations, which express 
the first order deviations due to the non-vanishing momentum transfer q: 

6(^S{A'^)) = S{-X) (59) 
6{Y{Ap),) = [S{X)J_ 
6[SFiP-A-\P/2-l))) = SFiP/2 + l)[Si-X),^/2-i^SFiP/2 + l), etc., 

we can consider the limit g ^ 0. After using the homogeneous Faddeev equation Eq. (|87D , 
we are left with the eight diagrams depicted in Fig.^. Note that the equality Fig.p.(a) -|- 
Fig.|.(6) ( = Fig.|.(d)) = Fig.| holds only in the limit g — 0, which is emphasized 
in Fig.^ by "g 0". By comparing both sides of Eq. (|^) in the limit q ^ 0, the 
explicit spin-parity projection]^ shows that the first two diagrams are proportional to 

^{a)a^p^_^M^^^(/3)fe('p^^ contribute to qa, and the other diagrams are proportional 

to Ii(°)"(P)-^75S(-A)m('^)''(P), which contribute to ^^^vtv- Note that the 1st, the 3rd,the 

4th and the 5th diagrams in Fig.^ come from the quark current contribution[Fig.^(a)], 
and the others come from the exchange current contribution[Fig.§.(b)]. 



5.2 Numerical Results in the Chiral Limit 

We first have to explain the choice of the parameters. There are four parameters: the 
cutoff A (Euclidean sharp cutoff), the current quark mass tjiq, the effective coupling 
constant in the pionic channel and the effective coupling constant in the qq scalar 
diquark channel Qsd- We fix the first three parameters (A, mo, (7,1) by solving the gap 
equation Eq.( |101| ) for the constituent quark mass M, and Eq.( |107| ) for the pion mass 771^^ 
and decay constant /,r, imposing the following three conditions: (1) m^r = 140 MeV, (2) 
= 93 MeV, (3) M = 400 MeV. The resulting values are A = 739 MeV, = 10.42 
GeV~^, mo = 8.99 MeV. Once these parameters are fixed, we consider the chiral limit by 
taking ttiq keeping A, g^, fixed. Eq.( [101| ) provides the mo dependences of M = M(mo). 
We treat gsd as a free parameter (independent of mo) to generate different nucleon masses. 
Note that, in the case m^^ = 140 MeV, gsd/ Q-k = 0.66 gives the experimental value of 
the nucleon mass 171^ = 940 MeV. For convenience, we plot the nucleon mass m^r (for 

^*The spin-parity projection is performed by using the helicity formalism. (See for detailH) 



17 



gsd/dn = 0.66), the quark- diquark threshold (mA? + mgd) where mgd is the scalar diquark 
mass, the pion mass m^, the pion decay constant and the constituent quark mass M 
against the current quark mass mo in Fig.|10|. We also plot the nucleon mass in the chiral 
limit against gsd/gn [solid line], together with the physical case(m^ = 140 MeV) [dashed 
line] in Fig.|Tl]. 

Our results for Qt^nn and in the chiral limit are plotted against gsd/dn by the 
solid line in Fig.|T2| and Fig.^, respectively. We obtain q^nn = 13.2 and qa = 1-32 for 
9sd/9ir = 0.66 compared to the experimental values g^^N = 13.5 and gA = 1.26. We 
will discuss the extension off the chiral limit in the next subsection. 

To estimate the violation of the GT relation, it is convenient to define a quantity: 

A f-Kg^NN /„„N 

Ag = . (60) 

mNgA 

In the chiral limit, the deviation of Ag from 1 is solely due to the cutoff artifact. We plot 
Ag against gsd/gn in Fig.|n| in the chiral limit [solid line]. It is seen that the violation of 



the GT relation is up to 4 %. In particular, for the reasonable case gsd/gn = 0.66, the 
violation is only 2 % in the chiral limit. 



5.3 PCAC Violation off the Chiral Limit 

In the chiral limit, Ag works as the measure of the violation of the PCAC relation due 
to the cutoff artifact. However, off the chiral limit, Ag does not work as the measure any 
more, since it contains, in addition to the unphysical cutoff artifact which we are going 
to estimate here, the physical effect of the non-vanishing current quark mass. Therefore, 
in order to estimate the cutoff artifact, we need to construct a quantity which picks up 
only the cutoff artifact for any values of rriQ > 0. To this end, we introduce another form 
factor m(q'^) through the following relation: 



#(a;)75y^(a;) 



Nl{L) = e^^X(^) ^75ir u^L) x ^^(g^). (61) 



Eq.(pTD generalizes in the following way: 

2mNgA{q^) + q^hAiq^) = 2moiA{q^). (62) 
In particular, in the limit 0, since now m^r 7^ 0, we have 

2mNgA = 2mou(g^ = 0). (63) 
Now we define a quantity Ap as follows: 

Ap = = (64) 

niNgA 
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By construction, the deviation of Ap from 1 is solely due to the cutoff artifact for any 
value of rriQ. To evaluate iAil"^), we first have to solve the following equation for the 
operator insertion on a constituent quark line similar to Eq.( 



5SP{x,y) 



T 



Sy{x,z) M75y 1 Sy{z,y) 



(65) 



+ J d'z'SP{x,z' 
I 

E2#Tr 



6SPiz',z') 
5p^{z) 



Sf{z\y). 



r r. 



Since (7ax = 0, we can parameterize 6T,/6p^ as 



6p^{z) 



«75- 



I{q^) then satisfies the following equation: 



with the solution 



2z^7^n55(g')-2^^7^n55(g')/(g'), 
-2z(7.n55(g') 



(66) 



(67) 



(6^ 



l + 2«^?,n55(g2)- 

Note that this quantity of course has the pion pole. The evaluation of iA{<f = 0) amounts 
to the 2nd to the 8th diagrams given in Fig.|^, where q^H{q^) is replaced by 1 + I{q^). 



In Fig.|T5|, we plot Ap against the current quark mass ttiq for gsd/ 9tt = 0.66 case [dashed 
line]. It is seen that the cutoff artifact is again within 4%. 



5.4 The Off-shell Qt^nn off the Chiral Limit 

We next evaluate gnNN off the chiral limit (m^ = 140 MeV). In this case, whereas the 
definition of gA = gA{.<f = 0) does not change]^, since the pion pole of hA^q^) is shifted, 
g-KNN is extracted according to: 



hA{q^ 



''^firg-K 



NN 



q2 _ ^2 



+ 



(69) 



In order to extract g^NN from the nucleon matrix element of the axial current, it is neces- 
sary to evaluate the form factor /i^(g^). However, at this stage it is still difficult to evaluate 
hA{(f) for non- vanishing momentum transfer in the relativistic Faddeev approach^. We 

in the case of 



15r 



^The evaluation of the iso-vector qa off the chiral limit was already done in ref. 
5ax = 0. 
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The numerical evaluation of the on-shell g^NN off the chiral limit is currently under consideration. 
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confine ourselves to evaluate the off-shell gnNNil"^ = 0) = 9nNN defined by the single-pole 
dominance approximation to hA^q^) as 

2 

~g^^^ = -fhAiq' = 0). (70) 

Note that, whereas the nearest cut in the physical fiAiq"^) is the three pion cut (g^ > 9m^), 
the Cutokosky rule(cf. p. 315 in |T5[) suggests that the nearest cut in our hA^q^) is the 
qq cut (g^ > 4M^). This is mainly due to the mean field approximation for the vacuum, 
and due to also the leak of the confinement in the NJL model. Since 4M^ is larger than 
9m^, the single-pole dominance approximation may work better in our case. We note 
that, although we could subtract the qq cut contributions from H{q^), it is impossible to 
subtract it from the remaining part, because the calculation refers only to g = 0. The 
explicit evaluation shows that the subtraction of the qq cut contributions from H{q^) leads 
to only a small difference. We note that the quantity Qt^nn is not only one of the possible 
off-shell extensions of the on-shell g^NN, but also the value of the axial form factor hA{q^) 
at g^ = up to the well-defined numerical factor given by Eg. ([7D|) . This enables us to 
derive an analytical expression of Ag which follows from the PCAC relation by neglecting 
the small cutoff artifact. (See Appendix 0.) 

We plot Qt^nn against gsd/ dn in Fig.|l^ [dashed line], and qa in Figjl^ [dashed line]. 
The value of cjnNN is 13.5 for Qsd/ 9^ = 0.66. This value is quite reasonable compared 
to the experimental value. We also plot Ag off the chiral limit [dotted line] and Ap 
[dashed line] against Qsa/ g-n in Fig.^. The reader might suspect why the validity of the 
GT relation could be improved by going off the chiral limit. The reason is that the effect 
of non- vanishing tjiq works into the opposite direction compared to the UV-cutoff artifact. 
To see this, we plot Ag against the current quark mass mo [solid hne] in Fig.|T3| for the 



case gsd/ g-K = 0.66. We also plot Ap against mo [dashed line], which is used to indicate 
the size of the cutoff artifact contained in Ag. It is seen that Ag is a monotonically 
decreasing function of mo. However, the experimental data Ag = 1.06 suggests that, as 
far as we believe that = 140 MeV is really close to the chiral limit, Ag should be 
an increasing function in the vicinity of mo — 0. We investigate this problem further in 
Appendix by taking into account also the effects of non-vanishing coupling constant 
in the isovector axial vector gg channel g^^x- (We leave this analysis for the appendix, 
because the non- vanishing g^.^ makes things quite complicated.) The main conclusions 
there are summarized as follows: (1) gA and cjnNN both decrease with increasing gs.x/gn, 
and increase with increasing mo. (2) The best fit of gA and g^NN could be obtained in the 
region < g^x/gn < 0.1, which, however, would depend on the quantity which we prefer to 
adjust. From this point of view, (7ax = is actually a rather good choice, because g-j^NN is 
very reasonable and gA = 1-33 is still close to the experimental value gA^^^ = 1-26. (3) For 
those values of g^^ which we examine, Ag remains to be a decreasing function of mo. (4) 
An analytic expression of Ag (Eq.( |137|) ) is derived by neglecting the small cut-off artifact 
and by assuming that the vacuum is approximated by the mean-field approximation. All 
the baryonic quantities disappear from this expression. In particular, this expression is 
valid even beyond the ladder truncation scheme for the 3qBS kernel. The discrepancy 
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between the analytic and the numeric Aq is due to the cutoff artifact. It is found to be 
within 3 %. (5) As a consequence, to make Aq to be an increasing function of mo and to 
obtain Aq — 1.06, we have to go beyond the vahdity of this analytic expression of Aq. 
Therefore, all we can do is either to improve the gap equation for the vacuum beyond the 
mean field approximation or to estimate the on-shell g^NN- We do not further investigate 
this problem in this paper. 

6 Summary and Discussions 

In this work we reviewed how to evaluate expectation values of quantum one-body oper- 
ators in the framework of the 3qBS /Faddeev equation by introducing classical external 
fields as a technical tool. In the 3qBS approach, the expectation values are obtained 
by sandwiching, the functional derivative of the 3qBS kernel with respect to the cor- 
responding external field between the 3qBS amplitudes. In the Faddeev approach, the 
expectation values are obtained by sandwiching, between the Faddeev amplitudes, the 
functional derivative of the Faddeev kernel with respect to the corresponding external 
field. We gave the criterion for 3qBS kernels to give rise to the PCAC relation correctly. 
For practical purpose, we also gave the sufficient condition for the validity of this criterion 
by introducing the local "axial" gauge transformation of the external fields. We applied 
the sufficient condition to several 3qBS kernels. The main results are as follows: (1) If 
the 3qBS kernel is truncated in the ladder truncation scheme keeping only the qq inter- 
action in the scalar diquark channel, the PCAC relation is obtained correctly. (2) If the 
qq interaction in the axial-vector diquark channel is included, it is necessary to include 
also the qq interaction in the vector diquark channel to give rise to the PCAC relation 
correctly. (3) Even if the qq interaction due to the qq exchange in both the pionic and the 
sigma mesonic channel is included, the correct PCAC relation is obtained. Concerning 
the point (2), we note that the vector diquark channel is often considered to be not im- 
portant. This is because the non-relativistic quark model suggests that the contribution 
from this channel to the nucleon mass is suppressed in the non-relativistic limit. However, 
to respect the chiral SU{2)l x SU{2)r symmetry, the qq interaction in the vector diquark 
channel should be included, even if it is expected to give a negligible contribution to the 
nucleon mass. We should note, however, that the relativistic Faddeev equation including 
all the qq interactions in the scalar, axial- vector and vector diquark channel amounts to a 
two-dimensional integral equation with 14 x 14 matrix structure even after the spin-parity 
projection is carried out, which requires a tremendous effort to be solved. 

Although these truncation schemes give rise to the PCAC relation correctly, the reg- 
ularization scheme, which has been adopted in practical numerical calculations so far, 
does not respect the chiral symmetry, leading to the violation of the Goldberger-Treiman 
relation. To estimate this violation, we carried out the numerical evaluations of Qa and 
g-^NN in the chiral limit in the simplest case, i.e., the ladder truncation scheme keeping 
only the qq interaction in the scalar diquark channel. We found that the violations are 
up to 4 %. In particular, for the case gs/g^K — 0.66, which reproduces the experimental 
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nucleon mass, the violation is only 2 %. The value of g-^NN in the chiral limit is 13.2 
which is quite close to the experimental value 13.5, and qa becomes 1.32 compared to the 
experimental value 1.26. We next estimated the PCAC violation due to the cutoff artifact 
off the chiral limit. We found that this violation is again within 4%. 

In the relativistic Faddeev method, it is still difficult to calculate form factors for non- 
vanishing momentum transfer, which is needed to extract the on-shell Qt^nn off the chiral 
limit. So we defined the off-shell Qt^nn by means of the single-pole dominance approxi- 
mation to hA^q^)- Although we obtained a very reasonable result g-irNN = 13.5 for the 
case Qsd/gn = 0.66, the effect of non- vanishing current quark mass mo on the Goldberger- 
Treiman violation {AQ{mo)) was found to be in the "wrong" direction: Whereas the 
experimental value of Ag is 1.06, which suggests that Aq should be an increasing func- 
tion of mo in the vicinity of mo = 0, our AG(mo) is a decreasing function of mo. We 
tried to resolve this problem (Appendix 0) by taking into account the effect of non- 
vanishing effective coupling constant in the iso-vector axial-vector mesonic channel g^x- 
The main results are as follows: (1) qa and g-rrNN both decrease with increasing g^x/gn, 
and increase with increasing mo. (2) The best fit of gA and g-rrNN could be obtained in 
the region < gax/ g-K < 0.1, which, however, would depend on the quantity which we 
prefer to adjust. From this point of view, (7ax = is actually a rather good choice, since 
g-KNN = 13.5 is a very reasonable result and gA = 1-33 is still close to the experimental 
value gA^^^ = 1-26. (3) For those values of g^x which we examined, Aq remains to be a 
decreasing function of mo. (4) An analytic expression of Aq was derived by neglecting the 
small cutoff artifact and by assuming that the vacuum is approximated by the mean-field 
method. All the baryonic quantities disappear from this expression. In particular, this 
expression is valid even beyond the ladder truncation scheme for the 3qBS kernel. The 
discrepancy between the analytic and the numeric Aq is due to the cutoff artifact. It was 
found to be within 3 %. (5) As a consequence, to make Aq to be an increasing function of 
mo, we have to take into account the effects which are beyond the validity of the analytic 
expression of Aq. Hence, all we can do are either to improve the vacuum beyond the 
meanfield approximation or to estimate the on-shell g^NN- 

We finally give a comment on the iso-scalar g^^\ It is straightforward to extend our 
formalism to the chiral U{1)l x U{1)r case. The Ua{^) anomaly in QCD is simulated in 
the NJL model as an explicit Ua{^) symmetry breaking. It is easy to see that, whereas the 
vector, axial- vector and tensor diquark channels form closed chiral U{1)l x t/(l)_R singlets 
separately, only a combination of the scalar diquark and the pseudo scalar diquark forms 
a closed chiral U{1)l x U{1)r doublet. To isolate the t/A(l) breaking contribution, one 
can parameterize the two coupling constants gsd and gpd as gsd = ^ + ^A, gpd = —A + 6X. 
Now, in the ladder truncation scheme, it is only 5X that can provide the Ua{^) breaking 
effects to 3qBS amplitudes, because A and the other qq interactions respect the t/A(l) 
symmetry. The qq interaction in the pseudo-scalar diquark channel is often considered 
to be irrelevant from the non-relativistic analogy. However, setting g^^ = corresponds 
to a particular choice of t/A(l) breaking, i.e., SX = gsdf^, and this particular choice is 
not based on any of the underlying physics of Ua{^) breaking. Due to this reason, we 
suggest that the qq interaction in the pseudoscalar diquark channel should be included 
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for a reasonable estimate of the iso-scalar g^^\ even if it is expected to give a negligible 
contribution to the nucleon mass. 
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Appendices 



A The Faddeev Equation 

The aim of this appendix is to summarize the notations of the relativistic Faddeev equation 
and to provide the derivations of some of the relevant relations involving the Faddeev 
method which are essential to the other parts of this paper in a self-contained manner. 
In this appendix, summations over repeated indices are not implied, and, unless explicitly 
indicated, the external fields are understood to be absent, i.e., v = a = 0, m = mo. 



A.l The Relativistic Faddeev Equation^ and the Green's Func- 
tion 

We begin with the 3qBS equation (see Eq.(^)) in the absence of the external fields: 

G = Go + KG; K = R\ + K2 + K3, (71) 

where the index i of Ki refers to the spectator quark. The formal solution of G is expressed 
by using the resolvent of K as follows: 

G = ^Go. (72) 

Note that the resolvent exists in a mathematical sense. However, because K is an un- 
bounded operator, it is difficult to interpret it as it stands. Therefore, we adopt the 
Faddeev prescription. We introduce the following Faddeev decomposition of the Green's 
function: 

G = Go + G^ + G^ + G^; G' = KiG. (73) 
We insert the following resolvent identity of K into Eq. ([7^) : 

We obtain 



1_A' 1-A', ')\- 



which is further inserted into the defining relation of G* in Eq.(^). We are left with the 
following closed equations for G^,G^,G^ (the Faddeev equations): 

^'-T^«« + r^SG^ (76) 



^^A rather good pedagogical introduction to the Faddeev equation is found in [ p7[ . 
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It is possible to simplify these coupled equations into a single closed equation for by 
using the identical particle nature of the three quarks. and G^ are obtained from G^ by 
means of simple permutation operations. For this purpose, it is convenient to introduce 
the cyclic permutation operator which is defined as follows: 

(Z?/')(xi,X2,X3) = V(x2,X3,Xi). (77) 

Note that Z is easily implemented by using delta functions. We give a list of obvious 
relations: 

^3 = 1, Z{1 + Z + Z'^) = {1 + Z + Z'^)Z = 1 + Z + Z^, (78) 
Ki = Z'K^Z-\ ZG = GZ = G, ZGq = GqZ = Gq, G' = Z'G\ 

Now Eq. (^) reduces to the following closed integral equation for G^ (the reduced Faddeev 
equation): 

= -r^Go + -^{Z + Z^)G\ (79) 
i - A3 i - -K3 

Kz 

which is depicted in Fig.|^.(a). The two-quark resolvent — is depicted in Fig.^(b). 

1 — A3 

(cf. Eq.(||)) It is not so hard to identify the so-called "Z-diagram" structure (the quark 
exchange diagram), which is provided by a combination of the permutation operator Z 

(or Z'^) and two external quark propagators of the two-quark resolvent — . We refer 

1 - A 3 

to the kernel of this integral equation as the Faddeev kernel Kp: 

K,^-^{Z + Z\ (80) 
1 - A3 

3 



Now the formal solution of G to Eq.(|79D is given as: 

By inserting this into the first relation in Eq. ([73|) , we have another formal representation 
of the Green's function G: 

G = G„ + (l + Z + Z=)^^G„. (82) 

A. 2 The 3qBS Amplitude and the Faddeev AmpUtude 

To obtain the form of the Green's function G near the nucleon pole, we first diagonalize 
the Faddeev kernel Kp, regarding the total four momentum as a parameteij^: 

KFi)[n]p\ = \n{p^)^[n;p\, ip[n]p\KF = \n{p^)il^[n]p\. (83) 



^The numerical procedure to solve this homogeneous Faddeev equation is explained into detail in 
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We adopt the following normalization condition: 

V^K;p>[n;p] =K(p2)(27r)V^)(p'-p)5„,.; = \ (84) 

where A^(p^) = d\n{p^) / d{p^) . Now the Faddeev kernel is expressed as: 

KF = Y.j (27r)Wn(p^) ^"^^^^^^^'-^^^^^'-^^' ^^^^ 

Note that, in reasonable cases, the eigen- modes associated with the largest eigenvalue are 
four-fold degenerate — they correspond to the nucleon (the ground states in the sector of 
the baryon number 1). The degeneracy is due to the iso-spin 1/2 and spin up/down. We 
refer to these eigen-modes as the (off-shell) eigenvalues and Faddeev amplitudes for the 
nucleon and denote them as \n{j)^), '>P[N^{p)] and iIj[N^{j))\. The nucleon mass is 
obtained by solving the following equation: 

Ajv(p' = ml) = 1. (86) 

The associated eigenvectors satisfy the homogeneous Faddeev equations: 

mm] = Y%<z + z')i;[Ni{L% m:{P)] = i^[K{P)]j%:{z + Z% (87) 

For simplicity, we suppress the time components of the total four momenta in the on-shell 
Faddeev amplitudes. We define the 3qBS amphtudes for nucleon states by: 

^[K{P)] = (1 + Z + Z^)ij[N^{P)] (88) 

In order to see that these definitions of the 3qBS amplitudes are reasonable and that the 
normalization scheme adopted in Eq.(|^) is consistent with the covariant normalization 
of the ket vectors in Eq.(^, we first consider the form of the Green's function near the 
nucleon pole as follows (cf. Eq.(|82D): 

^ - ^'^^^^^hJ (2vr)W.(P^) l-A.(P^) l-K,^' ^''^ 



E 



d'P ^ 1 ma{P)MNa{P)] 



(27r)W^(m^) -X'j,{m%) - + le 

I d^P niN 1 mN''(P)mN^(P)]. 

{27iY EN{P^)P,-EM{P^)+r^ i aK )\ [ aK )\ 
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This result is consistent with the expression which is expected from the canonical operator 
analysis. Note that the above manipulations are exact at the nucleon pole. 

Next, we show that these 3qBS amplitudes satisfy the homogeneous 3qBS equation 
Eq.(PD: 



^ Z''K3Z-\1 + Z + Z^)^[N^{P)] 

j=l,2,3 

(1 + Z + Z^)Ks{l + Z + Z^)i;[N^{P)] 



(90) 



1 + Z + Z' 



[z + z')ij[N:{p)] = ^KiP)] 



1 + Z + Z^ 



(91) 



:i+z+z^)K,- 



K. l + Z + Z^ ^ 



I-K3 3 

where, to obtain the 3rd lines of Eq.(0) and Eg . (|9l|) , we used the following identities: 

1 



'1 + Z + Z' 



(92) 



[1 + Z + Z')iIj 



■{Z + Z')^ + {Z + Z')^lj 



[Z + Z')^lj, 



which immediately follow from the homogeneous Faddeev equations Eq. (p7|) . 



A. 3 The Matrix Elements in terms of the Faddeev Amphtude 

Here, we derive an explicit expression for the matrix element of the axial vector current 
operator in terms of the Faddeev amplitudes, which is, in the practical applications, more 
convenient than the expression in terms of the 3qBS amplitudes. To avoid cumbersome 



notations, we introduce a shorthand notation: Oj.^ = 



Sa''{x = 0) 



(i=l,2,3) and 01 



O^.^, where the index i refers to the spectator quark. Note that O^.^ = Z^O^.^Z' 

j=l,2,3 

Now we have from Eq.(|T 



iV„»(P) A^ix = 0) iV^(L) 
= t[K(P)]Ojvl/[iV^(L)] 



maip)] 



l + Z + Z')0%^{l + Z + Z'mNl{L)] 
1 1 



1-7^3 



■(Z + Z>[iV^(L)] 
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where, to obtain the fourth hne, we used Eg. (p2D . Since the permutation operator Z does 
not depend on the external fields, 5/5a|^ only hits one the constituent quark propagators 

in the two-quark resolvent — in Kp. (cf. Eq.(|80|)) 

1 - K 3 

We comment on the ladder truncation scheme. In this case, the qq interaction is 
point-like. By combining the diagramatic expressions of K^^ [Fig.^.(a)] and the two-quark 

resolvent — [Fig.|^.(b)], we see that the diagrams involved in ibKpjba"' are classified 

1 - A 3 

into two types, i.e., (1) 5/5a'^^ hits one of the internal quark propagators (i.e., in the 
qq bubble diagram) in the ladder sum, (2) S/Sa'J^ hits one of the two external quark 
propagators. With the aid of this classification, a straight forward diagramatic argument 
shows that the diagramatic expression of Eg . (|9^) is given by the three diagrams Fig.|^.(a), 
Fig.||.(b) and Fig.||.(c). The first type leads to Fig.||.(c), which we refer to as the "diguark 
current" contribution, and the second type leads to Fig.||.(a) and Fig.H.(b), which we 
refer to as the "guark current" contribution and the "exchange current" contribution, 
respectively. 



B Notations of the NJL Model 

The aims of this section are to define some of the relevant guantities with explicit exam- 
ples, and to summarize the notations of the NJL model. In this section, the external fields 
are understood to be absent, i.e., v'^{x) = aJJ(x) = 0, m^x) = mo. 



B.l The Elementary Effective Coupling Constants and the Closed 
Chiral Multiplets 

Here we define the elementary effective coupling constants in qq and qq channels. We 
start with the Lagrangian density Eg.(|l]). For definiteness, we virtually distinguish the 
two V^'s and the two t/^'s from each other, respectively, i.e., 

>^:/ = E^?r(^irv>2)(^3r7/'4). (94) 

r 

The elementary interaction, which is depicted in Fig.|l^.(a), can be classified into the fol- 
lowing three types: (i) qq direct channel [Fig.|16|.(b)], (ii) qq exchange channel [Fig.|l6[(c)], 
(iii) qq diguark channel [Fig.|l6|.(d)]. They are related by Fierz identities to each other. 
We define the effective coupling constants (yi^^^dir)^ ^(ggexch) g(qq) ii^Yongh. the following 
relations: 

J^i = T.^2i29S'^{^^^c.r.M2){^3r^nM,) (95) 

a i=0 A=0 
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3 8 



3 8 



a i=OA=0 



Ta^T"' are Dirac gamma matrices, where a runs over S(scalar), V(vector), T(tensor), 
A(axial- vector) and P (pseudo-scalar), i.e., 

(rQ)a=s,v,T,A,P = (l^T/.^o-Mi', 7/^75, 75) (96) 

(r°)„=s,v,T,A,p = (l,7^^^^ 757^75). 

Ti {i = 0,1,2,3) arc the iso-spin Pauh matrices with tq = 1. They are normahzed 
according to tr(Tjrj) = 26ij. Note that, both in qq and qq representations, tq corresponds 
to the iso-scalar, and Ti {i = 1, 2, 3) correspond to the iso-vector channels. Pa are the 

R 

rescaled Gell-Mann color matrices, i.e., /3o = 1, Pa — y 2"^^ ~ 1)2, ■■■,8) with 

the normalization tr{PAPB) = ^Sab- In representations, Po corresponds to the Ic 
mesonic channels, and Pa foi A = 1,2, ■■■ ,8 corresponds to 8c mesonic channels. Note 
that Pa for A = 2, 5, 7 are anti-symmetric and the others are symmetric. Therefore, 
in qq representation. Pa for A — 2, 5, 7 correspond to 3c diquark channels, and Pa for 
A — 0,1,3,4,6,8 correspond to 6c diquark channels. We define the effective couphng 
constants g^*^^^ as follows: 

_ ^('jgdir) (gqexch) /q^x 
9aiA — 9aiA + 9aiA ■ ' ) 

We consider two examples (i) the original NJL type interaction Lagrangian: 

Li = g - E(^75TiV)') , (98) 

where Ti is the isospin Pauli matrix, and (ii) the color current type interaction Lagrangian: 

= E (v^TmY V') , (99) 

where Aq is the color Gell-Mann matrix. We give a list of these effective coupling constants 
giig) — g,(99<iir) _|_ giqqexch) g{qq) thcsc two interaction Lagrangians: 
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Original NJL Type 


S 


V 


T 


A 


P 


le (gg), 1 = 


^ + ^/12 


0-9/6 


+ 5/12 


0-5/6 


+ 5/12 


Ic (qq), I = 1 


0-9/12 


+ 


0-5/12 


+ 


-5-5/12 


8, {qq), 1 = 


+ 9/12 


0-5/6 


+ 5/12 


0-5/6 


+ 5/12 


8c (qq), I = 1 


- g/12 


+ 


0-5/12 


+ 


- 5/12 


3c (gg),/ = o 


9/Q 


-9/12 


9/Q (!) 


-5/12 (!) 


5/6 


3c (gg), / = 1 


0(!) 


9/12 (!) 





5/12 


0(!) 


6e (gg), / = 


9/Q (!) 


-5/12 (!) 


5/6 


-5/12 


V6 (!) 


6c (gg), / = 1 





5/12 


0(!) 


^/12 (!) 





Color Current Type 


S 


V 


T 


A 


P 


Ic (gg), 7 = 


+ 2^/9 


0-9/9 


+ 


+ 5/9 


- 25/9 


Ic (gg), / = 1 


+ 2^/9 


0-9/9 


+ 


+ 5/9 


- 25/9 


8c (gg), / = 


- 9/36 


-9 + 9/72 


+ 


- 5/72 


+ 5/36 


8c (gg), / = 1 


- 9/3Q 


+ 9/72 


+ 


- 5/72 


+ 5/36 


3c (gg),/ = 


9/Q 


-9/I8 


0(!) 


9/18 (!) 


-5/9 


3c (gg), 1 = 1 


9/9 (!) 


-9/I8 (!) 





5/I8 


-5/9 (!) 


6c (gg), / = 


-^/18 (!) 


5/36 (!) 





-5/36 


5/I8 (!) 


6c (gg), 1 = 1 


-9/18 


5/36 


0(!) 


-5/36 (!) 


5/I8 



For completeness, we also presented the effective couplings in the 6c diquark channels, 
which do not contribute directly to the color singlet baryon states. Note that the two 
■0 fields and the two t/j fields are originally undistinguished Grassmann fields. Therefore, 
in qq representation, unless Ta'y5C~^TiT2PA is anti-symmetric, the contribution vanishes. 
Due to this reason, half of the effective coupling constants in the qq channels actually 
vanish, and these cases are indicated by "(!)" in the list. 

Due to the chiral SU{2)l x SU{2)r symmetry of the interaction Lagrangian, it can 
be already seen from the above list that some of the effective couphng constants are 
grouped together. In fact, the straight forward apphcation of the chiral SU{2)l x SU{2)r 



30 



transformation leads us to the following list of "closed cliiral multiplets"|^: 





A/ T — n\ 
— oj 








A/ r — 

^^K, i — UJ 








(at r\\ 

[A,I = 0) 








/ /I r n\ 

/ = Oj 








(iSJ = 0) 


-(PJ = 


1)) 




(iSJ = 0) 


-(PJ = 






f(Q T - I) 


- (P I - 


0)) 




((Q T - l) 


- (p I - 


0)) 






— (T T — 






[{1,1 — V) 


— (T T — 






(iV I =1) 


~ (A I = 


1)1 




((V I = 1) 


- fv4 / = 


1)) 




(s, / = o) 








{s,i = i) 








(t,/ = i) 








(r,/ = o) 








(p,/ = o) 








(^,/ = i) 








((F,/ = 0) 


-(A,/ = 






((^,/ = i) 


-(A/ = 


0)) 



Note that the total number of independent coupling constants of the chiral SU{2)l x 
SU{2)Fi symmetric NJL type interaction Lagrangian is at most eight. Therefore, if we fix 
all the eight coupling constants in the qq channels, all the effective coupling constants in 
the qq channels are obtained as their linear combinations. In the ladder truncation scheme 
of the 3qBS kernel, we can fix the coupling constants in the 3c (qq) channels based on 
the calculations for color singlet baryons. However, the coupling constants in the 6c (qq) 
channels remain free. One can then use these remaining coupling constants (related to 
the qq coupling constants) to reproduce the mesonic properties. 

To avoid cumbersome notations, we introduce the following abbreviations: 

i/TT — ypj=i,i^ — y5,/=o,ie' "sd — ys,/=o,3c' "^x — yA,i=i,ic ~ yv,/=i,ic- \^^^) 

B.2 The Vacuum of the NJL model at the Mean Field Level 

Here, we summarize some of the notations of the NJL model in the vacuum and the 
mesonic sectors. 

B.2.1 5(ax = Case 

The gap equation is given by: 

^^For the precise meaning, see Section. 4. 
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which provides the spontaneous breaking of the chiral symmetry in the NJL model at 
the level of mean field approximation. M is the constituent quark mass. The pion mass 
and the pion decay constant are obtained from the two-point axial current correlator as 
follows: 







TAl{x)K{0) 







q2 _ ^2 _|_ 

Here we adopted the following definition of the pion decay constant: 



contmuum. 







■ipx 



(102) 



(103) 



where the covariant normalization (ti°'{p) T:^{k)j = 2ym'^+fP(2Tr)'^5'^^\'p — k)5ab is 
adopted. The explicit form of the axial current correlator can be easily obtained in 
the ladder approximation. However, in order to emphasize the consistency between the 
vacuum, the mesonic sector and the baryonic sector in our formulation, we prefer to use 
the external field method. Since the canonical operator expression of SSp/Sal is 



idiS^jx, z)ap 







T^„(a;)^^(z)4(y) 







(104) 



the explicit form of the axial current correlator is obtained in the following way(cf. eqs. 



/ 



^4^gi5(x-y) /q 



j d^xe^^^^-J'^Tr 



TAl{x)Al{y) 
( 

(7m75i^) 

V 



(105) 



SSp (x, x) 

Kiy) 



1^ r u 1 A n5A(9^)(-22^0H5A(g^) 

•T^VWJ^ab ~ -iq/j-qi^^ab . , r,- — rr r 2\ 

4 4 1 + 2«5(^H55(g^) 



where 



n55(g')^, 



ab 



qtM'n5Aiq'^)Sab = -J 



H^i/(g)5, 



ab 



(2^ 
d^k 

(2^ 
d^k 

(2^ 
d^k 



Ti({^,T-)iSp{k + q){ri,T'')iSF{k) 
Tr ((75r'^)z5p(A; + q){^^^^T')iSF{k) 
Tr ((757/.T^)^5i.(fc + g)(75r")z5p(A;) 



(106) 



Tr (7^75r")«S^(A; + q){l^lurySp{k) 



Jab 



(27r)4 

^1 Mq') + 5ab L.-^ 



n2(g'). 
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The explicit expressions for m^r and /^r are obtained from: 



(107) 



where U'^^{q ) 



B.2.2 ^ Case 

Even if we switch on g^^^ ^ 0, the gap equation does not change. However, the axial 
current correlator changes in the following way: 



(108) 



Jab 



g2 



{q'H{q^)U,A{q') + G{q^)n,{q') 



1 (q^,q,\N{q^) 

-dab ' ' 



where H[q^), Gi(g^), G2{q^) and D{q^) are defined in Appendix 0, and N{q^) is given as 
follows: 



N{q^ 



--^q\-2ig^) (n5A(g'))' (l + 2ig^Jlr{q^ 
-i(l + 2z^,n55(g')) {2ig,Jii{q^ 



(109) 



The explicit expressions for and f^^ are obtained from: 



= D(m2), 



Ar(m2 



(110) 



where D {q ) 



dD{q^ 



B.3 The Regularized Bubble Integrals 

The regularized expressions for the bubble integrals with the sharp Euclidean cut-off are 
as follows: 



^^^^ ' J {2tiY Jo {kl + - q^x{l - x)f 



[111) 
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n,,(,=) = 24«f g|/;..^j^^^^,-L-^^ (112) 

'^^^^ i (27r)4 7o (A;| + M2 - g2a;(i _ a;))2 V ; 
These expressions satisfy the following identities: 

n5A(o) = 2Mn'55(o) (114) 

q^n^M) = 2Mq^n,A{q^) (115) 

ni(o) = n2(o). (116) 

We comment here on the expression for n^j,(g). In order to obtain the above expression 
for Il^^{q), we have to use the following prescription: 



' dx \, -TT^ (117) 



(27r)4 Jo (M2 - P - g2a;(i _ ^^y 

This is due to the following reason. The straight forward application of the sharp Eu- 
clidean cut-off leads, rather than to Eg. ( P-15|) , to the following identity: 

q'U.^iq) = -q'n%\q) + 2Mq,U,Aiq'), (118) 

where Il^^\q) is 

^Z\Q)Sat = - / ^Tr ((7,r«)z5^(fc + g)(7.r^)^5^(fc)) . (119) 

Since the non-vanishing longitudinal part of n|^^)(g) is an unphysical cutoff artifact, which 
spoils the significance of the outputs, the prescription Eq.( |117| ) is often used to suppress 
it. (See p. 178 in ref. |18[.) Note that, if the dimensional regularization is applied. 



the expression Eq. (|117|) vanishes identically. Now the analogous subtraction should be 
performed on the l.h.s. in Eq.( |118| ), otherwise the Gell-Mann-Oakes-Renner relation would 
be terribly violated in the case (/ax 7^ 0. Note that, once we use this prescription, the Gell- 
Mann-Oakes-Renner violations are less than 1% for m^r < 140 MeV. 



C Matrix Element in terms of the 3qBS/Faddeev 
Amplitudes in the Rest Frame 

In principle, we can use any Lorentz frames to calculate the Lorentz invariant form factors 
of a matrix element. But in practice, it is convenient to make use of the rest frame, 
because this frame is often used in the calculation of the mass and the 3qBS amplitudes. 
The aim of this appendix is to give expressions for the matrix element in terms of the 
3qBS /Faddeev amplitudes in the rest frame by using the Lorentz covariance. Unless 
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the opposite is explicitly indicated, the external fields are understood to be absent, i.e., 
V = a = 0, m = niQ. To make the Lorentz transformation of the spinor simpler, we select 
the final momentum P = (mjv,0) and the initial momentum = {K~^)^^P^, where 
A'^jy = represents the boost. We introduce the Fourier transforms of the 3qBS 



amplitudes and of 



KM 



vl/[iV|(L)](M') {2iTf5'^'\L - L) 

= I dSidS2d'y^e-'^'^y^nmL)]iyi^y2,y^) 

' d^xid^x^d^x^ t[iV^(P)](xi,X2,X3) e^^P'"^ 



(120) 



d'^xid^X2d^X3 / d'^yid'^y2d'^y3 e 



0';,{P,L;p,p',l,l') 

d^P r d^L 



P, L are the total momenta and p,p',l,l' are the relative momenta, which are defined 
though the following relations: 



P = J2Pi P = m - (1 - v)iPi + P2), p' = v'pi - (1 - V')P2 
L = Ek i = vh-ii- V)ih + h), I' = v'h - (1 - v')k, 



where < 77,77' < 1 are arbitrary real numbers. We note that, unlike the nonrelativistic 
approaches, the choice of r] and 77' is almost completely arbitrary in the relativistic quan- 
tum field theory^". The delta functions in the first two relations in Eq. (|120|) are due to 
the translational invariance. It is important to note that the Jacobians associated with 
the variable change, i.e., (^1,^2,^3) ^ [P-iP-iP') and {h^h.h) ^ {L,l,l') are 1, i.e.. 



d^Pi d'^p2 d'^p^ d'^P d'^p d^p' 



d% d% d% 



(27r)4 (27r)4 {2tiY {2tiY (27r)4 (27r)4 ' (27r)4 {2tiY {2tiY 
Now the matrix element reduces to the following expression: 



d^L dH dH' 
(2^(2^(2^ 



K{p) 4(0) Ni{L)) nK{p)] o;. nm^)] 

d^p rfV f dH dH' 



(123) 



(27r)4 (2vr)4 y (27r)4 (27r)4 
X t[iV„" (P)] (p, p') OiiP, L- p, p\ 1,1')^ [m{L)] (/, /') 



•"For complete expositions, see pi 
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rfV f dH dH' 



(27r)4 (27r)4 y (27r)4 (27r)4 
X t[iV,"(P)](p,pO Ol{P,L-p,p'HH') S(A-^)^3vl;[iV^(P)](A/,Ar) 
d^p rfV /■ t^^^ dH' 



(27r)4 (27r)4 7 (27r)4 (27r)4 
X t[Ar^(P)](p,p') Oj(P,L;p,p',A-i/,A-i/') S(A-i)®3vI;[iV^(P)](/, /'), 

where the second equahty follows from the following Lorentz transformation property of 
the 3qBS amplitude of the rest frame: 





Ti,{xr)i,{x2)i^{x^)\Nl{K-^P)) (124) 
[T^(xi)^(x2)^(x3)] S{K~^)\nI{P) 
S(A"i)®2 (O \T'4){Kxi)%1){Kx2)'4){Kx^)\ N^{P) 

^ is the pure boost matrix, ^(A^^) is the boost operator in the 



the boost matrix in the Dirac 



where (A-^)^'' = (e"^ 

Fock space, S(A-i) = e^(-^), (5(-A) = ^(t^^(-A)'^^) 

bispinor space, and S{tC)®'^ = S{A) ® ^(A) ® S{A). In the last line of Eq.([T^, with the 
aid of the Lorentz invariance of the Jacobian, we rotate the integration variables /, /' by 
the boost A, i.e., / ^ A~^l, I' i— A^^l'. The problem thus reduces to the "matrix element" 
calculation of Oj(P, A^^P; p, p', A'^l, A-^l')S{A-^)^^ between two 3qBS amplitudes in the 
rest frame. 

By using the homogeneous Faddeev equation Eq. (|87D , the Faddeev amplitude is ex- 
pressed by the 3qBS amplitude as follows: 

V^[Ar„"(P)] = ir3^K(P)]. (125) 

Hence, the Lorentz transformation property of the 3qBS amplitude implies the following 
Lorentz transformation property of the Faddeev amplitude: 

V^[Ar-(A-ip)](xi,X2,X3) = ^(A-i)«3^[Ar«(P)](Axi,Ax2,Ax3). (126) 

By repeating almost the same arguments, we are left with the following expression of the 
matrix element in terms of the Faddeev amplitudes in the rest frame: 

[K{p)\a'>^{x = o)\ni{l)) 

d^p d^p' r dH dH' 



(127) 



(27r)4 (27r)4 y (27r)4 (27r)4 
x^[Ar,"(P)](p,pOO^^,(P,L;p,p',A-\A"Y)5(A-^)^3^[Ar|(p)](z,r), 

where the Fourier transforms 'i/'[A^|(-P)](p,p')? '^WaiP)]iP^P') defined from iIj[N^{P)]{xi, X2, X3) 
■iIj[N^{P)]{xi, X2, X3) similar to Eq. ( |12(]| ), and Op {P,L]p,p',l,l') is defined from: 



Op {xi,X2,Xs;yi,y2,y3) 



i6Kp{xi,X2,X3;yi,y2,y3) 



6al^{x 



0) 



(128) 
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D Non- vanishing ^ax 
D.l General Formalism 

The aim of this appendix is to consider how to extract qa and g-^NN in the case (^ax 7^ 
providing an analytical expression of the GT violation AG'(mo). Instead of repeating the 
argument similar to the one given in Section which would become quite lengthy, we 
derive the expression of the form factors in a different manner, which would be easier for 
the readers to understand. 

We first have to consider functional derivatives of the constituent quark propagator, 
i.e., SSp/Sa^^ and dSp/Sp^, which are obtained as the solutions to Eq. (|34|) and Eq. (|65D . 
We parameterize 612 /Sa^ and 6T,/6p^ as follows: 



(129) 



r 



= h{q' 



2 ^ 



'q^4 T^" 



G2{q )[Y r 7517 



n 



+ I2 



n 



H{q^)^ Gi(g^), G2{q^), Ii{,q^) and hiq"^) satisfy the following coupled equations: 



G2{q') 

H{q') \ 
Gi{q') I 

' hiq') \ 



2ig,Jl2{q^) + 2ig,Jl2{q')G2{q' 



(130) 



-2tg^U,A{q^) 
2tg^Jli{q'') 

-2ig^U,5{q^) 
2igaJJ5A{q'^) 



+ 



~2tg^U,,{q^) -2tg^U,A{q'') ^ 
q'^2ig^Jl5A{q'^) 2ig^Jli{q'^) j 

-2ig^Il,,{q^) q^{-2tg^)U,Aiq') 
2ig^Jl5A{q'^) 2ig^^Ui{q'^) 



H{q') 
Gi{q') 

( h{q') 



with the following solutions: 

-2ig^U5A{q'^) 



Hiq') 

Giiq') 
G2{q') 

hie) 



D{q' 



q^ {2ig^Jl,A{q^)) ( - 2ig^U,A{q^)) + (l + 2tg^U^,{q^)) {2ig,Jli{q' 



D{q') 

2ig^Ji2{q') 
I - 2ig^Ii2{q') 

1 - 22(7axni(g2)) - 2ig^Yl^^{q^)) + q" [- 2igJl^A{q^)) {2ig^J{^A{q^ 



Diq' 



2i^axn5A(g^ 

D(g2) 
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1 + 2zf7.n55(g') ) (1 - 2z(7axni(g2) ) - ( - 2ig^Ii^A{q^)) (2ig,Ji^A{q 



We define "baryonic parts" Bp{q'^) and B^iq"^) of tlie form factors through Fig.O. Now 
the matrix elements are expressed as follows: 



N{P) = 0)| N{L)) = u{pf- {gA{q^){l,l,) + hA{q'){q,l,)) u{L) 



(132) 



i^uiP)-^^^,u{L)j BA{q') (l + G^iq')) 
M l^-rm^^ (T^ R ( 2^ G'i(g2)-G2(g2) 



r 



+ luiP)j^ML)] Bpiq')q,Hiq' 



N{P) 



#(0)75y^(0) 



iV(L)^ = u{P)- {^Aiq'){^l5)) u{L) (133) 
(P)^Z75«(^) ) Bp{q^) (l + /i(g2)) + Lpfly^^^u^L)] BAiq')W), 



where the isospin and the helicity indices of the initial and the final nucleons are suppressed 
for simplicity. The form factors are expressed as 



BAiq') [l + G^iq' 
Bp{q-)H{q-) + 2m,BA{q')^l^^l^^^ 

q^ 

Bp{q^) (l + h{q^)) + 2mMBA{q^)h{q^). 



(134) 



9A{q') 
hA{q') 
^A{q') 

QttNN is extracted according to Eq. (|69D and qa by qa = dA^q'^ = 0). Note that, due to 
Eq.( |ll(j| ), Gi{q^) — G2{q^) is proportional to g^, which cancel q^ in the denominator in the 
expression of hA^q^)- 

By neglecting the small cut-off artifact, the PCAC relation Eg . (1281) leads to the rela- 
tion: 

2mNgA{q^) + q^hA{q^) = 2moM(g^), (135) 
which provides the following relation between Bp{q^) and BA{q^)'- 

2mr,BA{q^) (l + Gi{q^) - 2mo/2(g')) + Bp{q^) (q^H{q^) - 2mo(l + h{q% = 0. (136) 

By using Eq.( |131| ), Eq. (|134| ) and Eq. (|136| ), we are left with the following analytical ex- 
pression of AG(mo): 



AG(r«o) 



f-Kd-KNN 



rriNgA 

ml -2ig^U,AiO) 



(137) 



2iga 



2mo 1 - 2z(7axni (0) 1 - 2tgAll, (0) 



(n;(o)-n'2(o) 
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where ) = — , liolQ' ) = — . All the baryonic quantities disappear irom 

A(mo). Note that, once the mean-field approximation for the vacuum is adopted, this 
expression is valid even beyond the ladder truncation scheme of the 3qBS kernel, as far 
as the kernel satisfies our criterion of PCAC relation^. This, in particular, implies that, 
to improve AG(mo), we have to improve the vacuum or to evaluate the on-shell q-^nn- 



Therefore, the easiest way to improve the result of Ag(^«o) in Section is to improve 
the vacuum by including the effect of g^^ 7^ 0, since it is not so easy to treat the vacuum 
beyond the meanfield approximation. 

D.2 The Numerical Results 

We first have to explain the choice of the parameters. There are five parameters: the cutoff 
A (Euclidean sharp cutoff), the current quark mass mo, the effective coupling constant in 
the pionic channel g^^, the effective coupling constant in the iso- vector axial- vector mesonic 
channel g^.^ and the effective coupling constant in the qq scalar diquark channel gsd- We 
treat (7ax as a free parameter, and fix the first three parameters (A, mo, g-i^) by solving the 
gap equation Eq. (|101|) for the constituent quark mass M, and Eq. ( |110| ) for the pion mass 



m^r and the pion decay constant /^r imposing the following three conditions: (1) m^^ = 140 
MeV, (2) = 93 MeV, (3) M = 400 MeV. We fix ^^d by requiring that miv = 940 MeV. 
We give our numerical results for the three cases (1) gs.x/ g-K = 0, (2) ga^/ gn = 0.25, (3) 
gax/g-K = 0.5. The explicit values of the parameters are listed as follows: 

1- g^^Jgrr = 0: A = 739 MeV, mo = 8.99 MeV, g^ = 10.4 GeV'^, gjg^ = 0.66. 

2. g^Jgrr = 0.25: A = 812 MeV, mo = 6.99 MeV, g^ = 8.13 GeV'^ gjg^ = 0.694. 

3. g^Jg^ = 0.5: A = 874 MeV, mo = 5.75 MeV, g^ = 6.71 GeV^^ gjg^ = 0.72. 

Once these parameters are fixed, we consider the chiral limit by taking mo keeping 
g-iT, ga.x and gsd fixed. Eq.( p.01| ) provides the mo dependences of M = M(mo). We 
confine ourselves to non-negative values of gs.x/gn because of the following reasons: (1) 
our examples of the interaction Lagrangians give non- negative values, i.e., g^x/ g-w = 
in the original NJL type, and g^x/ g-w = 0.5 in the color current type. (2) As we shall 
see below, gA and g-j^NN increase with decreasing gax/g-K- If fi'ax < 0, we cannot adjust 
either g^ nor gnNN- (3) The increase of gA with decreasing gg,x is due to the fact that the 
negative g^x/g-K works as an attraction in the transversal iso- vector axial- vector mesonic 
channel (oi channel), which makes G2{q'^ = 0) to grow up, leading to the rapid increase 
of gA = gA{,<i^ = 0). (cf. Eq. ( |134| )) This "anti-screening" of gA does not seem to be 
reasonable [01. 



Now we present our numerical results. In Fig. 18, we plot g^NN against the current 



quark mass mo for these three cases of g^x- It is seen that gj^NN decreases with increasing 



^^Although the baryonic parts Bp{q^) and BA{q^) change, they stiU satisfy Eq.( 136 ). Hence Eq.(137) 
remains to be vahd. 
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Qax/g-n, and increases with increasing rriQ. The crosses in the figure are used to indicate 
mo which correspond to = 140 MeV. We use the diamonds to indicate mg which 
correspond to m^r = 2 x 140 MeV, which is used to indicate the vahdity of the single pole 
dominance approximation of Qt^nn- Note that the distance between and 4M^, i.e., the 
qq cut, is still 10 times larger than m^, as long as < 2 x 140 MeV. In Fig.|19|, we plot qa 
against mo for the three cases of (7ax- We see that qa decreases with increasing gi^/ g-K and 
increases with increasing mo. Note that, if ga.^/g-K = 0.1, then qa = 1-23 and g-^NN = 12.5. 
So the best fit of gT^NN and gA could be obtained in the region < g^^lgT, < 0.1, which, 
however, would depend on the quantity which we prefer to adjust. From this point of 
view, (j^ax = is actually a rather good choice. In this case, g^NN = 13.5 is very reasonable 
and gA = 1-33 is still close to the experimental value gA^"^^ = 1-26. Note that, due to the 
chiral symmetry, positive gax/gn implies an attractive interaction in the iso- vector vector 
mesonic channel (p meson). In order to describe p meson in the NJL model, however, we 
need stronger (jfax [@, 0, but, stronger ^fax leads to smaller qa and g-irNN in our calculation. 
(The situation could be improved by including qq interactions in the axial-vector diquark 
channel, which enhance g^NN In this case, one should also include the qq interaction 
in the vector diquark channel for the Goldberger-Treiman and the PCAC relation.) In 
Fig.^, we plot the numerical Aq together with the analytic Aq (Eq. ( |137] )) against mo. 
It is seen that they are monotonically decreasing functions of mo independent of ga^/g-K- 
We see that the discrepancy between the analytic and the numeric Aq is within 3 %. 
Recall that, when deriving Eq.( |137| ), we neglected the small cut-off artifact. Hence this 
discrepancy is solely due to the cut-off artifact. Note that the analytic Aq approaches 
1 as mo ^ in Fig.|2D|. In fact, based on Eq.CT), Eq.(|TT3D, Eq.(|TT§ and the Gell- 
Mann-Oakes-Renner relation, it is possible to prove that the analytic Aq approaches to 
1 as mo 0. To make Aq an increasing function of mo, it is necessary to go beyond the 
validity of the analytic Aq (Eq.( |137 | ) ). This implies that all we can do is either to improve 
the gap equation for the vacuum beyond the meanfield approximation or to evaluate the 
on-shell g^NN- In Fig. 21, we plot the PCAC violation Ap against mo for the three cases 
of gax/g-K- The deviation of Ap from 1 is solely due to the cutoff artifact. We see that the 
PCAC violations are only within 3%, which suggests the practical validity of our results. 



References 

[1] Y. Nambu, G. Jona-Lasinio, Phys. Rev. 122 (1960) 345; 124 (1961) 246. 

[2] U. Vogl, W. Weise, Prog. Part. Nucl. Phys. 27(1991)195. 

[3] S. Klevansky. Rev. Mod. Phys. 64 (1992) 649. 

T. Hatsuda, T. Kunihiro, Phys. Rep. 247 (1994) 221. 

R. Alkofer, H. Reinhardt, H. Weigel, Phys. Rep. 265 (1996) 139. 

C. Christov et al. Prog. Part. Nucl. Phys. 37 (1996) 1. 

[4] M. Kato, W. Bentz, K. Yazaki, K. Tanaka, Nucl. Phys. A551 (1993) 541. 
S. Kahana, G. Ripka, Nuch. Phys. A429 (1984) 462. 



40 



D. I. Diakonov, V. Yu. Petrov, P. V. Poboylitsa, Nucl. Phys. B306 (1988) 809. 
R. Alkofer, Nucl. Phys. B306 (1990) 310. 

Th. Meissner, F. Griimmer, K. Goeke, Phys. Lett. B 227 (1989)296. 

P. Sieber, Th. Meissner, F. Griimmer, K. Goeke, Nucl. Phys. A 547 459. 

[5] A. Buck, R. Alkofer, H. Reinhardt, Phys. Lett. B 286 (1992) 29. 

[6] N. Ishii, W. Bentz, K. Yazaki, Phys. Lett. B301 (1993) 165. 

[7] N. Ishii, W. Bentz, K. Yazaki, Phys. Lett. B318 (1993) 26. 

[8] N. Ishii, W. Bentz, K. Yazaki, Nucl. Phys. A587 (1995) 617. 

[9] H. Asami, N. Ishii, W. Bentz, K. Yazaki, Phys. Rev. C51 (1995) 3388. 

[10] M. Oettel, G. Hellstern, R. Alkofer, and H. Reinhardt, Phys. Rev. C58 (1998) 2549. 
H. Mineo, W. Bentz, K. Yazaki, Phys. Rev. C60 (1999) 065201 

[11] N. Ishii, Phys. Lett. B431 (1998) 1. 

[12] A. Bender, C. D. Roberts, L. v. Smekal, Phys. Lett. B380 (1996) 7. 
[13] M. Oettel, M. A. Pichowsky, L. v. Smekal, |nucl-th/9909082i 

[14] A. N. Kvinikhidze, B. Blankleider, Phys. Rev. C60 (1999) 044004, 044003, and 
personal communication. 

[15] C. Itzykson, J.-B. Zuber, "Quantum field theory", McGraw-Hill International Edi- 
tions (1985). 

[16] S. Pokorski, Section 12 in "Gauge Field Theories", Cambridge University Press 
(1990). 

[17] I. R. Afnan, A. W. Thomas, in "Modern Three-Hadron Physics" (Springer- Verlag, 
Berlin (1977), ed. A. W. Thomas). 

[18] K. Huang, "Quarks, Leptons and Gauge Fields". 

[19] M. Takizawa, K. Kubodera, F. Myhrer, Phys. Lett. B261 (1991) 221. 

[20] C. Fischer, Diploma thesis, Eberhard-Karls-Universitat Tiibingen (1999). 
R. Alkofer, personal communication. 



41 



Figure Captions 



The gauge transformation property of sP is depicted. A single line represents a 
propagator of the constituent quark in the presence of the external fields. The 
gauge transformation reduces to the multiplication of the phase factors Q{x) and 

niy). 



2. The functional derivative of the propagator, i.e., 



is depicted. A single 



line is a constituent quark propagator. Dashed lines are used to indicate the mo- 
mentum transfer q. Fig.|^.(a) corresponds to (?ax = case, and Fig.|^.(b) corresponds 
to (^ax 7^ case. 

3. (a) The 3qBS kernel in the ladder truncation scheme is depicted. A single line 
represents a quark propagator. A slash indicates an amputation, i.e., removement 
of the quark leg. 

(b) The diagram used to evaluate the nucleon matrix element for 3qBS kernel of the 
type Fig.^.(a) is depicted. "®" is understood in the sense of Fig.^. Diagrams of the 
same topologies are omitted. 

4. (a) The qq exchange improved 3qBS kernel is depicted. A single line represents a 
quark propagator. A slash indicates an amputation. A wavy line represents a ladder 
sum of qq bubbles as is depicted in the second line. Diagrams of the same topologies 
are omitted. 

(b) The diagrams with which to obtain the nucleon matrix element for 3qBS kernel 
of the type Fig.^.(a) are depicted. The second term is often referred to as the "meson 
exchange current" contribution. The diagrams of the same topologies are omitted. 

5. The proof that the ladder truncated 3qBS kernel gives rise to the PCAC relation 
correctly is diagrammatically explained. In the first step, the local gauge transforma- 
tion = e*'''^'^*-^^ is applied. Only the constituent quark propagators transform. 
In the next step, Eq.(H) is used. The last step is due to the fact that amputated 



propagators are delta functions. It is seen that the sufficient condition is satisfied. 

The proof that the qq exchange improved 3qBS kernel gives rise to the PCAC 
relation correctly is diagrammatically explained. In the first step, we applied the 
gauge transformation. In the second step, Eq.(|4T|) is used. Note that all the phase 
factors, which appear at the internal vertices, cancel themselves. The last step is 
due to the fact that the amputated propagators are delta functions. It is seen that 
the sufficient condition is satisfied. 



7. (a) The Faddeev equation Eq . (|79|) is depicted. A single line is a constituent quark 
propagator. A double line is a t-matrix in the qq diquark channel. (For precise 
meaning, see Fig.|^.(b).) (b) The diagramatic expression of the two-quark resolvent 
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— is depicted. A single line is a constituent quark propagator. A slash indicates 

1 — A 3 

an amputation. 

8. The diagrams which contribute to the matrix element {N\A^^{x)\N) are depicted. 
A blob with "N" which is followed by a triple line is a quark-diquark Faddeev 
amplitude of the nucleon. A single line is a constituent quark propagator, and a 
double line is a t-matrix in the qq scalar diquark channel. A dashed line indicates a 
momentum transfer g^. Fig.|.(a) is referred to as the "quark current" contribution, 
and Fig.H.(b) the "exchange current" contribution. Fig.|[(c) is referred to as the 
"diquark current" contribution. If the qq interaction is truncated to the scalar 
diquark channel, Fig.||.(c) does not contribute to qa, due to the iso-scalar nature 
of the scalar diquark. Fig.||.(a') is equivalent to Fig.|.(a), which is obtained by 
once iterating the Faddeev amplitude in the finial state by using the homogeneous 
Faddeev equation. By using Eq. (|127D and Y{^P)ii = S{A)^S{A-^), Fig.|.(a') and 
Fig.H.(b) reduce to the four diagrams depicted in Fig.||.(d). 

9. Diagrams which contribute to qa and Qt^nn are depicted. These eight diagrams are 
equivalent to four diagrams in Fig.H.(d) in the limit g — 0, which is indicated by 
"g —>■ 0" . The explicit spin-parity projection shows that the first two contribute to 
gA, and the others to QnNN- 

10. The nucleon mass rriN (for gsd/gpi = 0.66) (solid line), the quark-diquark threshold 
M + nisd {gsd/Qn = 0.66 case) (dotted line), the pion mass m-„ (dashed line), the pion 
decay constant (dot-dashed line), the constituent quark mass M (dot-dot-dashed 
line) are plotted against the current quark mass rriQ. The points which corresponds 
to = 140 MeV are indicated by the vertical dotted line. 

11. The nucleon mass is plotted against Qsa/ for the two cases (1) the chiral limit (solid 
line), (2) off the chiral limit [m^ = 140 MeV] (dashed line). The vertical dotted line 
indicates Qsa/ Q-k = 0.66, where rriN becomes 940 MeV for the case = 140 MeV. 

12. g-nNN is plotted against gsd/g-n for the two cases (1) the chiral limit (solid line), (2) 
off the chiral limit [171^^ = 140 MeV] (dashed line). The vertical dotted line indicates 
dsd/ diT = 0.66, where ttin- becomes 940 MeV for the case m^r = 140 MeV. 

13. The iso- vector gA is plotted against gsd/d-K for the two cases (1) the chiral limit (solid 
line), (2) off the chiral limit [m^ = 140 MeV] (dashed line). The vertical dotted line 
indicates g^d/ Q-k = 0.66, where mjq becomes 940 MeV for the case = 140 MeV. 

14. The violation of the GT relation Aq is plotted against gsd/d-K for the two cases (1) 
the chiral limit (solid line), (2) off the chiral limit [m^r = 140 MeV] (dotted line). 
The dashed line is the plot of the violation of the PCAC relation Ap off the chiral 
limit. The vertical dotted line indicates gsd/ g-w = 0.66, where rriN becomes 940 MeV 
for the case m^r = 140 MeV. 
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15. The violation of the GT relation Aq (solid line) and the violation of the PCAC 
relation Ap (dashed line) are plotted against the current quark mass ttiq for the 
case gsd/g-K = 0.66. The vertical dotted line indicates the current quark mass mo 
which corresponds to = 140 MeV. 

16. The interaction Lagrangian £/ is depicted in Fig.|TB[(a). The diagramatic interpre- 
tations of this interaction in the qq channel are classified into the two types: the 
direct channel [Fig.|TB[(b)] and the exchange channel (Fig.^.(c)). Fig.|r^.(b) is ob- 
tained from Fig.p!6|. (a) by using the Fierz identity. Another Fierz identity leads to 
the diagramatic representation of the interaction in the qq channel which is depicted 
in Fig.|16Kd). 

17. The baryonic parts Bp{q^) and B^iq'^) of the form factors are defined diagrammati- 
cally. A blob with "N" which is followed by a triple line is a quark-diquark Faddeev 
amplitude of the nucleon. A single line is the constituent quark propagator, and a 
double line is a t-matrix in the qq scalar diquark channel. A dashed line indicates 
a momentum transfer g^. 

18. g-KNN is plotted against the current quark mass mo for the three cases (1) g^x/ g-K = 
(solid line), (2) gaxlg-w = 0.25 (dotted line), (3) gs.x/ g-K = 0.5 (dashed line). The 
crosses indicate the points which correspond to m^^ = 140 MeV, and the diamonds 
indicate the points which correspond to m^^ = 280 MeV. 

19. gA is plotted against the current quark mass mo for the three cases (1) gs.^/ g-w = 
(solid line), (2) g^x/ g-K = 0.25 (dotted line), (3) ga.xl g-w = 0.5 (dashed line). The 
crosses indicate the points which correspond to m^^ = 140 MeV, and the diamonds 
indicate the points which correspond to m^ = 280 MeV. 

20. The numerical Aq is plotted against the current quark mass mo for the three cases 
(1) gax/gn = (sohd line), (2) g^^x/gn = 0.25 (dashed line), (3) g^x/gn = 0.5 (dot- 
dot-dashed line). The analytic Aq is also plotted for these three cases (1) by dotted 
line, (2) by dot-dashed line, (3) by dot-dash-dashed line. The crosses indicate the 
points which correspond to m^ = 140 MeV, and the diamonds indicate the points 
which correspond to 171^^ = 280 MeV. 

21. The violation of the PCAC relation Ap is plotted against the current quark mass 
mo for the three cases (1) g^x/ g-n = (solid line), (2) g^x/ g-n = 0.25 (dotted line), 
(3) ga.x/ g-K = 0.5 (dashed line). The crosses indicate the points which correspond to 
m^r = 140 MeV, and the diamonds indicate the points which correspond to = 280 
MeV. 
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Fig. 18 
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